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Abstract 

Nuclear anapole moments are parity-odd, time-reversal-even El moments of 
the electromagnetic current operator. Although the existence of this moment 
was recognized theoretically soon after the discovery of parity nonconserva- 
tion (PNC), its experimental isolation was achieved only recently, when a 
new level of precision was reached in a measurement of the hyperfine depen- 
dence of atomic PNC in ^^^Cs. An important anapole moment bound in 
^'^^Tl also exists. In this paper, we present the details of the first calculation 
of these anapole moments in the framework commonly used in other stud- 
ies of hadronic PNC, a meson exchange potential that includes long-range 
pion exchange and enough degrees of freedom to describe the five indepen- 
dent S — P amplitudes induced by short-range interactions. The resulting 
contributions of tt— , p—, and w— exchange to the single-nucleon anapole mo- 
ment, to parity admixtures in the nuclear ground state, and to PNC exchange 
currents are evaluated, using configuration-mixed shell-model wave functions. 
The experimental anapole moment constraints on the PNC meson-nucleon 
coupling constants are derived and compared with those from other tests of 
the hadronic weak interaction. While the bounds obtained from the anapole 
moment results are consistent with the broad "reasonable ranges" defined by 
theory, they are not in good agreement with the constraints from the other 



experiments. We explore possible explanations for the discrepancy and com- 
ment on the potential importance of new experiments. 



2 



I. INTRODUCTION 



The strangeness-conserving {AS = 0) weak nucleon-nucleon interaction is of consider- 
able interest. It provides the one experimentally accessible means of probing the neutral 
current component of the hadronic weak interaction, as this component plays no role in 
flavor- changing reactions. Furthermore the question of how long-range weak forces between 
nucleons are connected to the underlying elementary weak quark-boson couplings of the 
standard model is an important strong-interaction question, one with potential connections 
to poorly understood phenomena such as the A I = 1/2 rule. One of the challenges in 
the field has been the experimental determination of the various spin and isospin contri- 
butions to the low-energy weak A^A^ interaction, as this interaction is dwarfed by much 
larger strong and electromagnetic forces. The weak effects can be isolated only by precisely 
measuring tiny effects associated with the parity nonconservation (PNC) accompanying this 
interaction. Because the PNC effects are typically of relative size ~ 10~^, only one class of 
elementary NN scattering experiments, p + p, has reached the requisite sensitivity. PNC 
effects have also been isolated in nuclear experiments, but only a few nuclear systems are 
sufficiently well understood to permit theorists to relate the observable to the underlying 
A^A^ interaction. For these reasons there is interest in finding new experimental constraints. 

Shortly after Lee and Yang's proposal that weak interactions violate parity, Vaks and 
Zeldovich [0 noted independently that an elementary particle (as well as composite sys- 
tems like the nucleon or nucleus) could have a new electromagnetic moment, the "anapole 
moment", corresponding to a PNC coupling to a virtual photon. One contribition to the 
anapole moments of hadrons would thus arise from PNC loop corrections to the electromag- 
netic vertex. Despite some early work on the contribution of the nucleon anapole moment 
to high-energy electron- nucleon scattering 0, the interest in anapole moments might have 
been limited to theorists had not Flambaum, Khriplovich, and collaborators pointed out 
their enhanced effects in atomic PNC experiments in heavy atoms. As the anapole moment 
is spin dependent, it contributes to the small hyperfine dependence of atomic PNC. (The 
dominant PNC effects in such experiments arise from the coherent vector coupling of the 
exchanged Z° to the nucleus, and are thus independent of nuclear spin.) While nuclear-spin- 
dependent effects do arise from Vector(electron)- Axial (nucleus) exchange, this nuclear 
coupling does not grow systematically with the nucleon number A of the nucleus: naively, 
the axial coupling in an odd-A nucleus is to the unpaired valence nucleon. Flambaum et al. 
observed that the anapole moment of a heavy nucleus grows as A^/^, so that weak radia- 
tive corrections to spin-dependent atomic PNC associated with the anapole moment would 
typically dominate over the corresonding tree-level Z° exchange for sufficiently large A (A > 
20). This growth means that spin-dependent atomic PNC effects should be dominated by 
the anapole moment - a radiative "correction" - and measurable in heavy atoms. 

Nevertheless, spin-dependent atomic PNC effects are still exceedingly small, typically 
~ 1% of the size of nuclear-spin- independent atomic PNC effects. Despite considerable 
effort, only limits existed on the anapole contribution until very recently. However, with 
the Colorado group's measurement 0] of atomic PNC in ^^^Cs at the level of 0.35%, a 
definitive (Ta) nuclear-spin-dependent effect emerged from the hyperfine differences. This 
measuement is the principal motivation for the work presented here. The goal of the present 
study is to carry out an analysis of the ^^^Cs anapole moment that follows as closely as 
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possible the formalism developed and employed in other A^A^ and nuclear tests of the low- 
energy hadronic weak interaction That formalism is based on the finite-range PNC A^A^ 
potential of Desplanques, Donoghue, and Holstein (DDH), a potential that contains sufficient 
freedom to describe the long-range vr exchange and the short-range physics governing the 
five independent PNC S — P NN amplitudes The resulting vr— , p— , and tu— exchange 
PNC A^A^ potential is employed in estimating the loop contributions to the single-nucleon 
anapole moment and the exchange current and nuclear polarization contributions to the 
nuclear anapole moment for ^^^Cs. We also present results for Tl, where an interesting 
anapole limit exists [[7|,^. 

The current work extends the treatment of Ref. |^ by including heavy-meson PNC con- 
tributions, thereby going beyond long-range vr exchange to the full DDH potential. This 
extension is crucial in describing the isospin character of both the single-nucleon and nu- 
clear polarizability contributions to the anapole moment. The main results of our study were 
recently presented in a letter |TD[ . Here we give the technical details of the heavy- meson cur- 
rent and polarizability calculations, and discuss the associated shell-model calculations and 
their potential shortcomings. Our approach differs from most earlier calculations [^|lT]-[T5| by 
avoiding one-body reductions of the currents and potentials: exchange currents and polariz- 
abilities are evaluated from shell-model two-body densities matrices, modified by short-range 
correlation functions that mimic the effects of missing high-momentum components. We also 
use a form for the anapole operator in which components of the three-current constrained 
by current conservation are rewritten in terms of a commutator with the Hamiltonian, and 
thus explicitly removed. 

The paper is organized as follows. In Sec. II we define the anapole moment and the 
electron-nucleus interaction it induces, and discuss connections with the generalized Siegert's 
theorm. In Sec. HI we describe the DDH PNC A^A^ interaction arising from vr— ,p— , and 
u;— exchange and its connections with the S — P amplitudes. The treatment of the one-body, 
exchange- current, and polarization contributions to the anapole moment are given in Sec. 
IV. The summation over intermediate nuclear states in the polarizability is performed by 
closure, after calibrating this approach in a series of more complete shell-model calculations 
in lighter nuclei. Other technical details - particularly the rather complicated heavy-meson 
exchange- current evaluations - are presented in Appendices A through D. In Sec. V experi- 
mental values for the anapole moments of ^^^Cs and ^osxi are deduced from the corresponding 
hyperfine PNC measurements. Other tests of the low-energy PNC A^A^ interaction are dis- 
cussed, and the constraints they impose on various PNC meson-nucleon couplings described. 
We address the issue of uncertainties in the shell-model nuclear structure calculations and 
attempt to assess the effects of missing correlations phenomenologically. In the concluding 
Sec. VI we discuss the resulting discrepancies and possible future work that would help 
address some of the open questions. 

II. THE ANAPOLE OPERATOR AND CURRENT CONSERVATION 

In this section we describe the anapole moment in terms of a classical current distribution 
|]T^,[T^. The corresponding operator for a quantum mechanical current is obtained from a 
multipole expansion that satisfies the generalized Siegert's theorem. We illustrate, in a 
simple one-body nuclear model, the relationship between the anapole moment and the PNC 
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A^A^ interaction and the consequences of current conservation. 



A. Anapole Moments in Classical Electromagnetism 

Given classical charge and current distributions, p(x ') and j{x'), the scalar and vector 
potentials, and A{x), are obtain from integrals over the Green's function. After a 

Taylor expansion around the source point x one obtains 



AttIx — x'\ 



d'x 'p{x'){l -x' -V + hx' ■Vf + -- ■}-\-, (1) 

I 4:7T\X\ 



A{x)=ld'x' ^^"^ 



Att\x — x'\ 



d'x 'j{x'){l - f ' • V + ^(f ' ■ V)2 + ■ ■ -j-^. (2) 

In the scalar potential expansion, the first term inside the curly bracket generates the 
total charge (monopole) moment; the second term, the electric dipole moment; and the 
third term, a combination of the quadrupole and monopole charge moments fl^. For the 
vector potential, the first term vanishes as there is no net current. After carefully taking 
the constraints of current conservation and the boundedness of the current density into 
account (which place six constraints on the bilinear products j{x )iXj ), there remain three 
independent components in the second term, corresponding to the magnetic dipole moment 
of a classical current distribution. Similarly, the third term involves a symmetric product 
of two coordinates with the current, generating 18 independent trilinear combinations, with 
10 constraints. The remaining eight independent components comprise the static magnetic 
quadrupole moment and the El moment known as the "anapole moment" (AM). 

One can extract the vector potential due to the AM explicitly, 

1(™'^)(^) = i-a^ + —a ■ —)^--, (3) 

where 

a = [ d^x'x' X {x' X j{x')). (4) 
6 J 

(We multiply and divide by for consistency with the definition of a we will later in- 
troduce via the Dirac equation.) We can remove the second term in Eq. (^ by a gauge 
transformation, so that 

Current conservation allows Eq. (H) to be rewritten as 
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FIG. 1. A toroidal current winding generates a nonzero anapole moment. 



(6) 



M2 /■ 3 , ,2^.^,. 

a = — — J a X X j[x ). 

(We use the Lorentz-Heaviside unit in which a = e'^/ATrhc = 1/137.) Eq. (|^) is often 
presented as the definition of the AM P,p]2|- p!^ , p!9|j . However it is important to note that this 
form is obtained only after exploiting the constraints of current conservation. 

It is apparent, for the ordinary electromagnetic current, that the associated AM operator 
is odd under a parity transformation. Therefore a nonzero AM requires either the introduc- 
tion of an axial vector component into the current or a parity admixture in the ground state 
(allowing the ordinary electromagnetic current to have a nonvanishing expectation value). 
This requirement of PNC associates the AM with the weak interaction. 

Another important property is the contact nature of the AM vector potential. Thus 
an atomic electron interacts with the AM of the nucleus only to the extent that its wave 
function penetrates the nucleus. 

Figure |I| gives a classical picture of the anapole moment as a current winding. Although 
the currents on the inner and outer sides of the torus oppose one another, there is a net 
contribution because of the weighting (in spherical coordinates) of the current in the defi- 
nition of the AM, leading to an AM that points upward. The illustrated current distribution 
is odd under a parity reversal, as we have noted it must be for the ordinary electromagnetic 
current. If, however, the current has a chirality - a small "pitch" corresponding to a left- or 
right-handed winding that would signal PNC - a parity-even contribution to the operator 
would be induced. 



B. The Anapole Operator 

Although one could quantize Eq. (|^) directly to generate the anapole moment operator, a 
better procedure is to avoid the assumption of current conservation, as this is often violated 
in nuclear models. Switching to a standard spherical multipole decomposition yields the 
momentum-space charge and current operators [|20| 

PiQ) = E(-0'47rr;^(fi,)Mj-'(g); (7) 

J,M 
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TABLE I. Properties of multipole moments under parity and time reversal. A slash (no slash) 
denotes odd (even) behavior. 



T 



el 



rpmag 



fT 



PT 
PT 



M 



coul 



J = 

J = 1 

J = 2 
J = 3 



PT 
PT 

n 



J,M 

and the associated charge, transverse electric, and transverse magnetic multipole projections 
of definite angular momentum and (in the absence of PNC) parity 



MTmIq) = J d'x jj{qx)YjM{n^)p{x); (9) 

TUq) = I d'x -V X [jj{qx)Yfj,{Q^)] .Jix); (10) 
J q 

TZ'iQ) = J d'x Jjiqx)Yfj,in^) . Jix), (11) 



where q is the (outgoing) three-momentum transfer, jj the spherical Bessel function, Yjm 
and Yfji the ordinary and vector spherical harmonics, and T'^^(— 0g, —6q, the rotation 
matrix. 

The transformation properties of the possible multipole moments under parity (P) and 
time-reversal (T) are listed in Table |. Systems that are parity and time-reversal invariant 
can have only even-rank Coulomb moments (charge, charge quadrupole, etc.) and odd-rank 
transverse magnetic moments (magnetic dipole, magnetic octupole, etc.). The P- and T-odd 
moments, which would arise in the standard model from small CP-violating contributions 
to the weak interaction, correspond to the odd-rank Coulomb and even-rank transverse 
magnetic multipoles (electric dipole, magnetic quadrupole, etc.). The PNC but T-even 
moments, which would arise from the usual weak interaction, correspond to the odd-rank 
transverse electric multipoles, with the lowest of these being the dipole moment known as 
the anapole moment. 

For consistency with Eq. (^, we require 

V'A{x) = -j{x), (12) 
which then defines the anapole operator 



aiA=lim^^ Tf,. (13) 



7 



The simplest case is the general expression for the matrix element of a conserved four- 
current for a free spin-| particle 

- - ^^^^^?.75)f/(p). (14) 

from which the four moments of Table | can be immediately identified. The two vector 
terms define the charge Fi{q^) and magnetic -^2(0'^) form factors. The axial terms that 
follow are the anapole and electric dipole terms, respectively. The anapole term reduces in 
the nonrelativistic limit to 

iAq"^ - q 1^)15 ^ -^q \(T-qq-a) 



a{q 



2-\ 



M2 



q^a^, (15) 



showing that the current is transverse and spin-dependent. From this current we then have 
the AM operator for a nonrelativistic point particle 

a-ix = a(0)o-u. (16) 



C. Current Conservation and the Extended Siegert's Theorem 

The anapole operator aix has been defined in terms of Tf|, and it is well known that 
this operator can be transformed into other forms by exploiting the continuity equation. 
These forms are equivalent in calculations where consistent charge and current operators 
can be constructed and exact matrix elements evaluated. However we are interested in 
nuclear calculations where, when one goes beyond the simplest descriptions to models that 
treat the interactions among the nucleons, current conservation is not preserved. We lack 
a prescription for constructing the many-body currents consistently that are necessary for 
current conservation, and for addressing the renormalizations that account for the limited 
Hilbert spaces employed in nuclear models. In such cases there is a preferred form for Tf{, 
the form in which all components of the three-current constrained by current conservation 
are reexpressed in terms of the charge operator. 

A familiar example is the case of El transitions generated by the ordinary electromag- 
netic current. Then Tf| generates a one-body operator proportional to p/M, which is of 
order v/c, where v is the nucleon velocity. It can be shown that the exchange current 
contribution to Tf| is also of this order. As the exchange currents, in general, cannot be 
constructed faithfully, it follows that errors will arise that are necessarily of leading order in 
the velocity. 



Siegert pT| showed that the situation could be greatly improved by exploiting the con- 



tinuity equation 

V-j{x) = -t[H,p{x)] (17) 
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to write Tf, in the long- wavelength limit, entirely in terms of the charge operator. This 
generates the familiar dipole form of the transverse electric operator, proportional to tur, 
where u is the energy transfer. The importance of this rewriting is that the charge operator, 
which is of order (v/c)^, has exchange current corrections only of order (u/c)^, or of relative 
size ~ 1%. Thus the Siegert's form of the El operator is a far more controlled operator in 
nuclear calculations. 

A form of Tf| consistent with Siegert's theorem is in common use P2 



q J 

-^(^^f" j d'x jj^,{q^)Yfj^,,{n^)-M, (18) 

where = means the equality holds after taking matrix elements {f\0\i). This form has the 
correct leading-order behavior for transitions due to the first term, with the second term 
vanishing as g — >■ 0. But for a static moment, the first term vanishes; the leading-order 
behavior is then governed by the second term, which the naive Siegert's theorem does not 
properly constrain. 

However the extension of Siegert's theorem to arbitrary q was derived by Friar and 



Fallieros |]23|,p2[ : at every order in q those components of the current constrained by current 



conservation are identified and rewritten in terms of the charge operator. The result is 



J 



J + 2 J (2J+1)!! 
where gj and hj are polynomials in q |2^. Combining Eq. (^) and Eq. (|T^) one finds 



ax = -— J d^x x^[jx{x) + V2^{Y2{nx) ® jix))x] (20) 

This is the AM operator form used here and in our earlier work; to our knowledge all other 
analyses have been based on the naive form of Eq. (|^). 

We stress that the three transverse electric operators T^' ,T^' , and T^' are equivalent for 
simple one-body models which ignore nucleon-nucleon interactions, provided the resulting 
one-body currents are properly generated by minimal substitution. The differences in these 
operators arise when they are used in more realistic calculations. 



D. Simple Examples 

In this section we illustrate how this equivalence is manifested in non-interacting shell 
model calculations of the nuclear AM of ^^^Cs. The PNC interaction is also taken to be a 
one-body effective potential, Hp\;Q. 

The elements of the calculation include: 
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1. Extreme single-particle forms for the ground-state nuclear wave function. As ^'^'^Cs is 
an odd-even nucleus with J = 7/2+, the odd proton is placed in the lgT/2 shell, outside 
an otherwise fully spin-paired closed core. 

2. The strong Hamiltonian is a one-body harmonic oscillator potential with spin-orbit 
interaction. While this description is primitive, it does yield the proper ground-state 
spin and parity for the (nearly spherical) nucleus ^^^Cs. The harmonic oscillator wave 
functions allow analytic calculations of polarizabilities, etc. 

3. Hp\iQ is treated perturbatively: only linear terms are retained. 
Thus the resulting Hamiltonian is 

H^Ho + H^lc, (21) 

with 

^0 = E §P + lMu;'x{t f - fs{t) ■ m- (22) 

i=l 

Hpkc = E 2M ' ^ ' 

where uj is related to the harmonic oscillator size parameter h hy uj — l/Mb"^, the spin-orbit 
strength / can be determined from shell splittings near the 3s2dlg shell, and gs and gy, the 
isoscalar and isovcctor strengths in the one-body PNC potential, can be chosen to represent 
the average potential exerted by the core nucleons. The analytic expressions we obtain 
illustrate the functional dependence on all of these parameters. Thus we are not concerned 
here with specific numerical values. 
By minimal substitution 

H^H + e^ ; p^p-eA, (24) 
one can derive the charge and current densities to order 1/M, 

P(^) = ef:i±^5(3)(f-x-), (25) 



and 



1=1 



jconv{x) = eJ2—^^^^{Pi^^^''^\^ - ^i)}sym] (26a) 

Jma,{x) =eYj^^^^p^V X {a{t)6^{x - X,))- (26b) 

l.oX^) = ejz^-^^y{i) X <y{i)5^'\x - X,)- (26c) 
i=i ^ ^ 

Jp-vc(f) =ei:i±|i^«^.-Wi'='(S-f.). (26d) 

i=l 
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where the subscripts conv,mag, s.o., and PNC denote the current densities arising from 
convection (kinetic energy), magnetization (intrinsic nucleon spin), the spin-orbit interac- 
tion, and the PNC potential, respectively. The first three are vector currents while the last 
is axial vector. Current conservation is then easily verified 

V ■ (jconvix) + jmag{x) + js.o{x) + jpNc{x)) = -i[Ho + H^plfC , p{^)]- (27) 

Contributions to the AM are generated by the axial-vector current acting between the 
unperturbed nuclear ground state and by vector currents that contribute because HplfQ 
perturbs the ground state, 

^ Y^^ {HTf'^^\Xo){Xo\H?Nc\'^o) \ He) (28) 
XO -^V'o ~ -^xo 

where ipQ and Xos are single particle unperturbed eigenfunctions of definite (and opposite) 
parities, and the superscripts (A) and (V) label the components of Tf generated by the 
axial-vector and vector currents, respectively. 

The special case of no spin-orbit interaction is interesting because the first-order per- 
turbed wave function (in fact, the result can be generalized to all orders) is given by the 
Michel transformation [0] 

tpQ{x) i){x) = (1 - ig'^ ■ x)i)Q{x) 

^ e~--{9-)-^Q(f)^ (29) 

where g = gs ± gv for a proton(+) or neutron(-). Eq. (^) shows that Hp^Kj generates a 
spin rotation along the radial direction characterized by a small angle proportional to g and 
to the distance to the center of the nucleus. Consider an S1/2 state aligned along the +z 
axis. The spin probability around a ring, centered at the origin, would be uniform and in the 
+z direction: we visualize this as a uniform array of up-spinors. When the weak interaction 



is turned on, the Michel rotation will produce a spin helix [|1^ structure for this chain of 
spinors as shown in Fig. If we picture each spin as a small current loop, the combination 
of all horizontal spin components S*// can be viewed as a toroidal current winding producing 



an AM, as discussed in Sec. Ill A 



Moreover, if the Michel transformed wave function is used in a calculation of the AM, 
one finds that the contributions from jpNc and jconv cancel exactly, so that jmag is entirely 
responsible for the AM. Even with the inclusion of the spin-orbit interaction, the magneti- 
zation current remains the major contribution to the AM 

The sum over intermediate states in Eq. ( PBD simplies considerably in the harmonic 
oscillator since the momentum operator p only generates transitions of one huj. Thus the 
transitions that must be consider in the extreme single-particle limit are the simple Ip and 
2plh single-shell transitions of Fig. ^ [jl5[. 

Some preliminary algebraic manipulations are helpful. Using the commutation relation 

+iMu;V,f] = -z^, (30) 



'2M 2 ' ' M' 
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FIG. 2. Spin helix structure due to the parity mixing. 




Valence Excitation Core Excitation 

FIG. 3. Valence and core excitations produced by the PNC potential acting on an extreme 
single-particle harmonic oscillation ground state. The levels illustrated are appropriate for -^^^Cs. 



12 



the PNC one-body potential can be rewritten as 



H 



(1) 

PNC 



9 
2M 
■ 9 



a ■ p 



= «^([^o, a-x] + -[a-l,a-x]). 



(31) 



Using this result in the polarization sum yields 



2 



f J2{t\a ■ {x X r)\n){n\Tf''^\t)/{E, - E^) + H.c. 



(32) 



As a typical value for the nuclear spin-orbit strength is f /uj = a ~ 0.1, one can work to first 
order in /, yielding for the various AM contributions 

gM ^ 



a 



(conv) 



-{i\x a\i) 
o 

^{i\a ■ (x X r)\n){n\l{x'p- ix)\i)/{Ei - K) + H.c. , 

4 „ Z 



Jmag) 



gM 



{(z|(a-f) 



X — X a\i) 



a 



^{s.o.) ^ 



(PNC) 



gM^ 
8 

gM_ 
8 



f{i\x {{a ■ x)x — X (?)|i), 
{i\x'^a\i) . 



(33) 

(34) 

(35) 
(36) 



As the 0{p) terms from a^'^'") and a^^^^^ exactly cancel, a^"*"^^ determines the leading- 
order (LO) contribution 

gM 



-li{i\{a ■ x)x — X a\i) 



(37) 



For the next-to-leading-order (NLO) O(f^) contributions, wc approximate Ei — E^ — —hu 
(as the spin-orbit correction to this are higher order) and invoke closure 



(conv) 9-^ / - l f -> ^ -> ^, i\ ^ f 2-* -"MIA 

~ -a{i\\a ■ [x X l).,-[x p — x)\\i)., 



(38) 
(39) 



Therefore, assuming these matrix elements are of the same order of magnitude, one obtains 
for the relative sizes 



I (conv) I I 
\'^NLO IC'LO\ 

I (mag) i i i i 

Wnlo /(^Lo\ ~ |a| , 



a 



2/1 



{x^ 

62 



^1/3 



a 



2/1 



~ \a\ 



(40) 

(41) 
(42) 
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TABLE II. Matrix elements corresponding to four choices of the anapole operator. Note the 
dimensionless parameter a = f /u! = fMb'^, where b is the oscillator parameter. The results are in 
units of C = AnJ^^gs + gy)Mb'^ , another dimensionless quantity. 



a 



Tf 



-4/3 
l-4a 



Tf 



Cartesian 



conv 
mag 
s.o. 
PNC 
Tot. 



-1.1 I -5.4 
l-4a l+4« 
11.16 _|_ 50.22 



l-4a 



l+4a 



-15.6a I -30.6a 
l+4a 



l-4a 



13/2 



11.16-20O _|_ 50.22-9a 



l-4a 



l+4a 



l-4a 
11.16 
l-4a 
-8.8a 



l-4a 



11/2 



same 



-2.7 
l+4a 
50.22 
l+4a 
-19.9a 



l+4a 



+ 



-6 
l+4a 



11.16 _|_ 50.22 



l-4a 
-44a/3 



l-4a 



+ 



l+4a 
-33a 
l-|-4a 



22/3 



same 



0.5 

l-4a 



+ 



-10.125 
l+4a 



11.16 _|_ 50.22 



l-4a 



l+4a 



22a I -49.5a 
l+4a 



l-4a 



77/8 



same 



where in the last line we assume an odd-proton nucleus with A ~ 100, similar to Cs. 

In Table |I| we present our results for the AM of ^^^Cs in this single-particle scheme using 
the four different Tfs discussed previously (Tf,Tf ,Tf , and the one from the Cartesian 
decomposition using Eq. (|13D and Eq. (||)) to define our anapole operator. Agreement is 
achieved only when: i) all the cnvTents-conv , mag , s.o., and PA^C-and ii) a complete set of 
excitations- valence and core-are considered. This illustrates a point made earlier, that the 
use of incomplete current operators or Hilbert spaces breaking current conservation will in 
general lead to difficulties. 

The table also shows that the contribution of the magnetization current, which is sep- 
arately conserved (V ■ j^™"^) = 0), is independent of the choice of the anapole operator. 
This term is entirely responsible for the leading 0(a°) result (given the conv-PNC cancel- 
lation in this order). It is also apparent that the NLO contribution attributed to a given 
current depends on the anapole operator chosen: it is the sum over all contributions, not in- 
dividual contributions, that is kept constant in calculations satisfying current conservation. 
The numerical value of the LO contribution (61.4) is reduced by ~ 20% to 50.5 when the 
NLO contributions are included {a = 0.1), consistent with our earlier assertion that these 
corrections are perturbations. 



III. THE PNC NUCLEON-NUCLEON POTENTIAL 

The AM calculations presented here are the first to employ an A^A^ weak potential 
sufficiently general to describe long-range vr exchange and all five short-range S — P NN 
amplitudes. This section summarizes the isospin structure of the AS" = hadronic weak 
interaction and its description in terms of vr— , p—, and cj— exchange. 



A. The Isospin Structure of the Hadronic Weak Interaction 



The standard model specifies the weak charged and neutral currents, Jw and Jz, asso- 



ciated with the absorption/emission of weak bosons by quarks 
light quarks (u,d,s) are 



The couphngs to the 
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= COS 9cu-f''{l - -i^)d + sin Ocu-i^'il - 75)5 ; (43) 

1 8 - 4 

Jz = -7k ^{^T^'Cl - o sin' + 75)^^ - d-t^il - - sin' Ow + 75)^^ 

V^cosfc'^ o o 

-s7''(l-^sin2^H^ + 75)s}, (44) 

where Oc is the Cabbibo angle, with sin 6*0 ~ -22 and 6w is the Weinberg angle, with 
sin' 9w ~ 0.23. The effective quark-quark weak interaction at low energies can be described 
by a phenomenological current-current Lagrangian 

^Weak — -^{JwJw+JwJw'^JzJz), (45) 

By assigning proper isospin and strangeness quantum number to each quark field, we 
can decompose these hadronic currents 

Jw = cos + sin Ocj'^'^^] (46) 

Jz = Jr^ + Jr\ (47) 

where the first superscript denotes the change in isospin (A/), and the second in strangeness 

(A^*). The current J^'^'* drives the u ^ d transition, while J^k ^^ drives the u ^ s transition. 
We can construct the strangeness conserving (AjS^O) hadronic weak interaction Lagrangian 
density 

A^S=0) f , 2n T(l,0)t t(1,0) , . 2n T(|.l)t t(^,1) , rj \ 



An important aspect of this Lagrangian density is its isospin content. The symmetric 
product of two J]y'^^ currents forms A/ = (isoscalar) and A/ = 2 (isotensor) interactions, 

(-,0) 

while the symmetric product of two J^' currents forms a A J = 1 (isovector) interaction. 
Therefore the charged current weak NN interaction in the A/ = 1 channel is suppressed by 
tan' 9c relative to the A7 = or 2 contributions. As there is no isovector suppression for 
the neutral current, one concludes that the A7 = 1 A^'A'" channel provides experimentahsts 
their best opportunity for studying the neutral current component of the hadronic weak 
interaction. 

The physical states are strongly interacting composites, nucleons and mesons. The strong 
interaction dresses the underlying quark-boson couplings, and we have not yet developed the 
theoretical tools needed to evaluate the strong effects quantitatively. The physical couplings 
associated with the effective operators for nucleons and mesons are thus expected to differ - 
perhaps substantially - from the underlying bare couplings. One famous example of this is 
the AJ =1/2 rule in strangeness- changing weak decays: in experiments one finds a strong 
enhancement of AJ = 1/2 over AJ = 3/2 amplitudes, relative to expectations based on the 
underlying standard-model couplings and efforts to evaluate strong renormalizations. One 
reason for the interest in PNC is the hope that we can learn more about such strong effects 
by adding precise data on AS" = weak hadronic interactions. 
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B. Meson Exchange and the Long-range PNC NN Potential 

The most straightforward contribution to the PNC nuclear potential is from the direct 
exchange of H^^and between nucleons. Because of the small Compton wavelengths of 
these bosons (~ 0.002 fm) direct exchanges effectively occur only when two nucleons overlap. 
We do not yet have an adequate understanding of such short-range contributions to either 
the PNC or parity-conserving (PC) NN interactions. Fortunately, for energies characteristic 
of bound nucleons, the NN interaction takes place primarily at distances large compared to 
the nucleon size. This is due in part to the strong repulsion in the NN interaction at short 
distances, and in part because nuclei are moderately dilute Fermi systems. Thus we expect 
long-range contributions, which can be described without explicit reference to the structure 
of the nucleon, to dominate the PNC interaction at low energies. 

The strong PC NN interaction at low energies (< 400 MeV) has been quite successfully 
modeled in terms of meson-exchange potentials. The complicated short- distance quark and 
gluon dynamics governing this interaction are parameterized by various meson-nucleon cou- 
plings and phenomeno logical form factors constrained by experiment. This meson-exchange 
strong interaction model can be enlarged to include the weak PNC NN interaction by re- 
placing one of the strong meson-nucleon couplings by a weak coupling. All of the physics 
of W and Z exchange between quarks - and the attendant strong interaction dressing - 
is buried inside the weak meson-nucleon vertices. As in the case of the strong A^A^ inter- 
action, the weak vertices depend on momentum-independent meson-nucleon couplings and 
phenomenological form factors. For this model to make sense, one should, at a minimum, 
be able to derive a consistent and reliable set of meson-nucleon couplings from PNC observ- 
ables. Should such a set emerge, the longer term goal would be to develop a first-principles 
understanding of the relationship between the effective hadronic couplings and the under- 
lying standard model bare couplings, dressed by a complicated soup of strong quark-quark 
interactions. 

In developing a sensible meson-exchange model for the PNC A^A^ force, one must first 
truncate the tower of possible dynamical mesons, effectively "integrating out" those which 
do not contribute explicitly to the interaction. At small center-of-mass energies light mesons 
dominate the PNC potential because they have longer ranges. Candidates below the chi- 
ral symmetry breaking scale of ~ 1 GeV include the pseudoscalar mesons 7r(140MeV), 
r/(549MeV), and r/'(958MeV); the scalar mesons S(975 MeV) and 5(983MeV); and the vector 
mesons p(769MeV), u;(783MeV) and 0(lO2OMeV). One could also consider crossed two-pion 
exchanges, etc. Barton's theorem which states that CP invariance forbids any coupling 
between neutral J=0 mesons and on-shell nucleons, helps to restrict the possibilities, elim- 
inating exchanges of 7r°, 77, r/ , S, and 6^ (to the extent that CP violation can be ignored). 
Furthermore McKellar and Pick have argued that 5^ exchange can be regarded as a form 
factor correction to n"^ exchange [^, and is strongly suppressed relative to p and u. This 
motivates a PNC potential based on n"^, p°, p^, and exchanges. (We will present below 
another argument that will make this potential seem less arbitrary.) 

The PC and PNC meson-nucleon interaction Lagrangian density in the vr— , p— , and 
a;— exchange model is 

^pc = ig-KNNN 75f ■ vf A^ - QpnnN (7^ - 'i-^o-f.^q")? ■ p^N 
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-9.nnN\^^ - t^a^^qniO^N, (49) 

+]v'(/iV + ^L^3^^)7M75iV, (50) 

where gnNN, dpNN, and g^jMN are the strong tt— , p— , and a;— nucleon couphng constants and 
h^p'^''^\ and /i^^'^-* (the superscripts denote the rank of isospin) are the weak vr— , p— , and 
tu— nucleon couphng constants. (In the hterature f-^ is also frequently called or h^.) Note 
that the 75 convention is that of Bjorken and Drell, and that q is the outgoing momentum 
of the produced meson. (Both of these conventions are opposite in sign to those of 0.) 
Evaluating the one-boson exchange diagrams, where one of the vertices is PC and the other 
PNC, and making a non-relativistic reduction, one obtains the PNC NN potential 

H?lcir) = ^[r(l) X r(2)]3[<?(l) + <?(2)] ■ u^rl 

{{^'^^^^ ■ ^^^^ + yf^^^)^ + ^^^^'^ + ^^^^'^^^^'^'^^ - ^^^^ ■ 

X {(1 + fivWil) X a(2)] ■ u,{f) + P(l) - a(2)] . v,if)} 

+ [Go + ^[t-3(1) + r3(2)]} X {(1 + psWil) X a{2)] " ^..(^1 

+ [a(l) - a{2)] . Um + - T-3(2)][a(l) + a(2)] ■ [G^^ - F,t;,(f)]j, (51) 

where r = ri — r2, u{r) = [p, e~"^^ /Anr], v{r) = {p, e~"^^ / Airr} , and p = pi —p2- The various 
coefficients in this potential are products of PC and PNC couplings: = gnNNfn/V^, 
Fo = —gpNNh^p/"^, Fi = F2 = —gpNNhy2, Gq = —gwNNh'^/'^, and Gi = 

—gijNNh'l/2. We use the strong couplings g-^-NN = 13.45, gpNN = 2.790, and g^NN = 8.37. 
Vector dominance fixes the strong scalar and vector magnetic moments, ps = —0.12 and 
Pv = 3.70. Note that the vr-exchange channel is 1=1; numerically it dominates the isovector 
A^A^ weak interaction. This is the channel which tests the strength of the neutral current 
component of the hadronic weak interaction. 

While the field has seen considerable experimental progress in constraining the PNC 
meson-nucleon couplings, the theoretical situation has hardly advanced beyond the bench- 
mark analysis of Desplanques, Donoghue, and Holstein (DDH) , carried out twenty years 
ago. Using SU(6)iy symmetry, current algebra, and the constituent quark model, DDH 
related charged current components of /^r and the hy to experimental PNC amplitudes for 
AS* = 1 nonleptonic hyperon decays. Portions of the neutral current contributions were 
also related to hyperon decays, while the remaining pieces - unaccessible through symmetry 
techniques - were computed using explicit quark model calculations. Uncertainties associ- 
ated with the latter imply considerable lattitude in the theoretical predictions. The resulting 
"best values" and "reasonable ranges" are given in Table |T|. The case of /vr is particularly 
acute, as this coupling is nominally dominated by neutral current interactions. 

Subsequent to the DDH work, other approaches, such as soliton models |5DI] and QCD 



sum rules |^T[, have been applied to the weak meson-nucleon couplings. None of these 



approaches, however, has yielded a sharper theoretical picture. Part of the difficulty may 
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TABLE III. Weak meson-nucleon coupling "best values" and "reasonable ranges" (in units 
of 10"'') from the standard model calculations of Desplanques, Donoghue, and Holstein. For 
comparison, the last two columns give the corresponding results of Dubovik and Zenkin (DZ) and 
Feldman, Crawford, Dubach, and Holstein (FCDH). 



Coupling 



"Reasonable Range" 



"Best Value" 



DZ 



FCDH m 



hi 



0.0^11.4 
-30.8^11.4 
-0.38^0.0 
-11.0^-7.6 
-10.3^5.7 
-1.9^-0.8 



4.6 
-11.4 
-0.19 
-9.5 
-1.9 
-1.1 



1.1 

-8.4 
0.4 
-6.8 
-3.8 
-2.3 



2.7 
-3.8 
-0.4 
-6.8 
-4.9 
-2.3 



lie in the assumption of valence quark dominance for the hadronic weak interaction. In 
particular, it has recently been shown, in the context of chiral perturbation theory, that chiral 
corrections to the leading-order PNC nNN interaction may be large . These corrections, 
which have no analog in constituent quark models, reflect the presence of "disconnected" 
light qq sea contributions. Given the present interest of hadron structure physicists in the sea 
quark structure of light hadrons, the possibility of important sea quark contributions makes 

a particularly interesting object of study. Achieving agreement among all determinations 
of this coupling is, thus, important. As we observe below, the current interpretation of the 
Cs and Tl AMs in terms of DDH couplings shows that such agreement is not yet in hand. 

In can be argued that an analysis in terms of meson-exchange PNC couplings is in fact 
quite general, if limited to low-energy observables: the DDH couplings are a shorthand 
for another representation of the low energy PNC A^A^ interaction, one based on the five 
independent S — P amplitudes. The DDH description in terms of vr, p, and u exchange can 
be viewed as an effective theory, valid at momentum scales much below the inverse range of 
the vector mesons. At low momentum the detailed short-range behavior of the potential is 
not resolvable: thus one could characterize the vector-meson contribution to the weak A^A^ 
interaction by five strengths describing the five S — P amplitudes. A sixth parameter would 
be needed to describe vr exchange, as this interaction is long ranged. The six DDH couplings 
thus are equivalent to such a description of the weak potential. 

In an ideal world one would determine the low-energy A^A^ S — P amplitudes, or equiva- 
lently the six weak meson-nucleon couplings, by a series of A^A^ scattering experiments. Such 
experiments require measurements of asymmetries ~ lO"*^, the natural scale for the ratio 
of weak and strong amplitudes, inGprnl/gl^^. As we will detail later, only a single A^A^ 
measurement, the longitudinal analyzing power for Al for p + p, has produced a definitive 
result. This result has been supplemented by PNC measurements in few-body nuclei and in 
some special nuclear systems where nuclear structure uncertainties can be largely circum- 
vented, allowing the experiments to be interpretted reliably. An analysis of these results, 
which have been in hand for some time, suggests that the isoscalar PNC interaction - which 
is dominated by p and u exchange - is comparable to or slightly larger than the DDH "best 
value," while the isovector interaction - dominated by vr exchange - is significantly weaker 
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FIG. 4. One-body axial-vector currents contributing to nucleonic anapole moments are gener- 
ated by pion loop diagrams and by vector meson dominance diagrams. 
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FIG. 5. Two-body axial-vector currents. One meson-nucleon vertex is strong, the second is 
weak. 

0. As the isovector channel is expected to be enhanced by neutral currents, there is great 
interest in confirming this result. One reason for the interest in the ^^^Cs AM is the hope 
that spin-dependent atomic PNC measurements can provide such a cross check. 



IV. CONTRIBUTIONS TO NUCLEAR ANAPOLE MOMENTS 

The DDH meson exchange model - which we have argued provides a very general de- 
scription of the PNC A^A^ interaction at low energies - has become the standard formalism 
for discussing low-energy properties of the weak A^A^ interaction. We now extend this for- 
malism to nuclear AMs, discussing the various PNC meson-exchange mechanisms by which 
a virtual El photon can be absorbed by the nucleus. 



1. Fig. ^ illustrates a PNC pion cloud dressing of a nucleon (one pion-nucleon coupling 
is PNC and one PC) and a vector-meson pole graph, leading to El photon absorption 
by a nucleon. The axial currents corresponding to such pion-loop and vector meson 
dominance diagrams generate nucleonic AMs P,|33|-p^ , which we discuss in more detail 
in Sec. 



IV A 



The two-body Hp\j(j also generates two-body axial vector exchange currents (See Fig. 
^). The diagrams we evaluate include: i) pair currents, where El photons couple to 
the A^A" pairs excited by the two-body potential, and ii) transition currents, where El 
photons couple to the exchanged mesons. Detailed calculations are described in Sec. 
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FIG. 6. An opposite-parity polarization of the nuclear ground state induced by the PNC weak 
A^A^ interaction. 



3. The two-body Hpj^^ polarizes the nucleus, producing an opposite-parity ground state 
component. This component then couples back to the unperturbed ground state via 
the amplitude for absorbing a virtual El photon. The resulting polarizability requires 
one to sum over a complete set of opposite-parity intermediate states (Fig. |^). This 
is discussed in Sec. |IV Q . 

The dependence of these contributions on nucleon number A is important. As the one- 
body anapole contribution involves a coupling to spin, it is easy to see that the nucleonic 
contribution acts very much like a nuclear magnetic moment: in a naive picture of an odd- A 
nucleus as an unpaired nucleon outside of a spin-paired core, the core contribution cancels, 
leaving only the valence nucleon contribution. While that contribution will depend on the 
quantum labels of the valence orbital, there is no general growth of the nucleonic contribution 
with A. In contrast, it was the important observation that that polarization contribution 
grows as A^^^ ^ that led atomic experimentalists to realize that AMs might be measureable. 
This growth not only leads to larger AMs in heavy nuclei, but guarantees that the AM will 
dominate over other sources of spin-dependent PNC, such as direct V(electron)-A(nucleus) 
Z° exchange (another nucleonic coupling that effectively sees only the unpaired valence 
spin). Similarly, it was shown in ||^ that the exchange current contribution also grows like 
A^^^. Note that the polarization contribution could be additionally enhanced if the ground 
state is a member a fortuitous parity doublet. There has been some discussion of anapole 
(and electric dipole) moment enhancements because of such accidental near-degeneracies 

In Figs. - P the AM is shown interacting with an external photon. Yet the illustrated 
processes are not physical, as the anapole coupling vanishes for on-shell photons. The 
underlying physical processes involve a scattering particle - e.g., an atomic electron, the 
source of the virtual photon. It follows that the AM need not be a gauge invariant quantity: 
instead it is one of a larger class of weak radiative corrections - corrections naively of 
0{GF<y) - that together form a gauge invariant physical amplitude. Included in this larger 
set of radiative corrections would be various "box" diagrams corresponding to simultaneous 
exchange between the electron and nucleus of a photon and Z^, etc. However the long- 
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distance contributions to the AM of a nucleus - the meson cloud contributions and many- 
body contributions due to wave function polarization and exchange currents discussed here 



- are both dominant numerically and separately gauge invariant |]33| . This is one reason the 
set of anapole contributions associated with Hp^(j discussed here are of such interest. 

The calculations require wave functions for the nuclear ground state and one- and two- 
body transition density matrices for evaluating the effects of one- and two-body operators 
on the ground state. The wave functions were derived from shell model (SM) diagonal- 
izations with harmonic oscillator Slater determinants and with suitable residual two-body 
interactions. For ^^^Cs, the oscillator parameter is 6 = 2.27f and the canonical SM space 
is between the magic shells of 50 and 82, i.e., l(y'7/2-2(i5/2-l/iii/2-2(i3/2-3si/2- Calculations 
were performed with the five valence protons restricted to the first two of these shells and 
four neutron holes to the last three. This produced an m-scheme basis of about 200,000. 
Two interactions were employed, the Baldridge-Vary potential and a recent potential 
developed by the Strasbourg group both of which are based on the addition of mul- 
tipole terms to g-matrix interactions and are designed for the ^^^Sn region. As the results 
are very similar, here we only quote results from the Baldridge-Vary calculation. For ^''^Tl, 
an oscillator parameter h = 2.54f was chosen. The ground state was described as a proton 
hole in the orbits immediately below the Z=82 closed shell, i.e., 3si/2-2d3/2-2d^/2 (though 
the 1 /ill/2 lies between two d shells, we omitted this opposite-parity shell to keep the SM 
space manageable), and the two neutron holes are in the space between magic shells of 126 
and 82, i.e., 3pi/2-2/5/2-3p3/2-l«i3/2-2/7/2-l/i9/2- A simple Serber- Yukawa force was used as 
the residual interaction. 



A. Nucleonic Anapole Moments 

As illustrated in Fig. ^, the one-body PNC electromagnetic currents (parity even) can 
be derived from pion loop diagrams, where one meson-nucleon vertex is weak and PNC and 
the other strong and PC, and from vector meson dominance. After plugging these one-body 
PNC currents into Eq. (pO|), the one-body anapole operator takes the form 



A 

-body 



J2[as{0) + a,{0)TM^xi^)■ (52) 

i=l 

This form makes it clear that the contributions of spin-paired core nucleons cancel, leaving 
only the valence nucleon AM. The results from 0, where only the pion contribution was 
considered, are 

a,(0) ^ -0.193/^6, (53) 
a„(0) ^ -0.048/^e. (54) 

Thus the pion loops generate an isoscalar coupling that is about four times larger than the 
isovector one. Later this calculation was extended to include the p°-pole contribution by 



vector meson dominance 13311 . This work was further extended to included the full set of one- 



loop contributions involving the DDH vector meson PNC couplings [35], using the framework 



of heavy baryon chiral perturbation theory (HB^PT) and retaining contributions through 
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0(1/A^), where = AirFj^ ~ IGeV is the chiral symmetry breaking scale and = 93MeV 
is the pion decay constant. This yielded the nucleonic AM couplings 



as{0) 



-0.24/^ - 0.37hl - 0.11/1° + Omhl - lA3h\ + 0.0051/i„s-i^ 

71" 



+0M7{hf^ + ^^^^) - 0.3{hf + h'f) + 0.009Vi^ - 0.125/1^"^ , (55) 



^ -0.37(/i° + -^^^) - 0.12/ii + Omhi - lA3h\ + 0.9(/i°, + h^y) + 0.0051/i„2-ic 
e v6 4 



+0.047(-/if + -^^^) - 0.3(/i^'^ - Kf) + O.OOQVi^ - 0.125/if (56) 
v2 

The HBxPT result for the pionic contribution is consistent with the earlier pion loop esti- 
mates: the isoscalar coupling is 1.3 times the pion loop value, while the isoscalar coupling 
is zero to this order in xPT. However, the vector mesons greatly enhance the isovector AM. 
An evaluation using DDH best values shows that a^(0) ~ 7as(0). That is, the inclusion of 
the vector mesons enhances the AM and qualitatively changes its isospin character, with the 
proton and neutron AMs opposite in sign. The HB^PT calculation included non- Yukawa 
type irNN couplings (defined as h^s and h\s in [^) associated with derivative interactions. 
Here we include only the standard DDH contributions, omitting the rest. Using "best val- 



ues" for the neglected terms [^, this omission is estimated to generate a 3% error in the 
dominant isovector coupling and 100% in 0^(0). The reason for the omission is consistency: 
such derivative couplings are absent in the DDH PNC A^A^ potential, the parameters of 
which are constrained by experiment. A consistent treatment of the derivative coupling 
would require not only their propagation through the polarization and exchange current 
calculations for the AM, but also redoing the DDH potential fits to all other low-energy A^A^ 
and nuclear PNC observables. We leave this ambitious task to future work. 

Folding these expressions with our SM matrix elements {{I\\Eti^i^\\I) = -2-372 and 
2.532, {I\\Eti ^(0^(0 11^) = -2-305 and 2.282, for Cs and Tl, respectively) yields the results 
in Table IVl 



B. Exchange Currents 

The virtual El photon can also be absorbed on a pair of nucleons coupled by the PNC 
potential. Such PNC exchange currents are evaluated in the standard way. The transition 
matrix is derived and reduced nonrelativistically, retaining terms through 1/M. This re- 
sulting momentum-space current is then Fourier transformed to produce a coordinate-space 
two nucleon current. 




xf{k,pl -pi,P2 -P2), (57) 
where x is the field point, xi and X2 the source points. 
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In Appendix A we give the two-body charge and current operators in momentum 
space. In Appendix B we give the nonvanishing three-current coordinate-space operators to 
0(1/M), the forms needed for the AM calculation. The vr contribution, which turns out to 
dominate numerically, is 

3^'^\x, Xi, X2) = 3\Z.p^c\x, Xi, X2) + JjJp^Nci^, Xi,X2) + /('''''^^(X, Xi,X2) 

-eg^NNfn ..^^^ _ .^2) - r(l)3r(2)3)(a(l)5(=^)(f - xi) + a(2)5(3)(x - X2) 



Sv^vrM 

. Vi - a{2) ■ V2)\{6^^\x - xi) + 6^^\x - X2)){xi - X2) 



2 

^i^^^^ (3) _^ (3) le 

2 J J Xi — Xo 



'm^\xi—X2\ 
\Xi — X2I 



Even with the complete exchange currents in hand, evaluating their shell model matrix 
elements is a formidable task. (The one previous AM exchange current calculation treated 
only TT exchange .) For example, the form of j'^'''^ is far more involved than any of the pionic 
contributions. The procedure we follow is to first identify which currents are numerically 
significant by averaging the currents over the nuclear core. Once identified, full two-body 
evaluations are then performed for these cases. 



The one-body average, first performed for PNC potentials by Michel [^, involves direct 
and exchange terms 

(a|0«|/?) = E(«7|0(^)|/37) - («7|0(^)|7/?), (59) 

7 

where the sum extends over all single-particle core states. The averages are done in a 
Fermi gas, a simple choice because spin, isospin, and spatial averages can be performed 
independently. The nucleus is viewed as a single particle outside a spin-paired (but isospin 
asymmetric) Fermi sea. The one-body average operators are obtained in closed form, though 
the average done over the spatial functions produces, in general, a complicated but smooth 
function of the single-particle initial and final momenta (the Y and W functions below). 
The smoothness allows us to replace this function with an average value, with little loss of 
accuracy. Appendix C contains an example of this averaging procedure, while the full results 
for the various currents are listed in Appendix D. In the case of vr exchange the result is 

_ -pair) -; ^ , -^(tt -pair) ^ . :^(7r7r7)/-j -t \ 

J \x,Xi)—j^_pQ yx, Xi) -\- j^^piyc yx, Xi) -\- J yx,Xt) 

eg.NNU ^^^^ ^ ^ _ ^^^^^^ ^ p((iy'W)a5(3)(f _ f^) 



2V2Mml 
2 



ml 



iX^)pla5^'^ (f - X,) - (Fi) [a ■ V, [V, 5^'^ {x -x,)]] 



-\{Y2){a ■ V„ {V„ 6^'\x-x,)}}] (60) 

The one-body estimate of the exchange current contributions to the AM can be obtained 
by plugging the averaged currents into Eq. (^). The Fermi-gas-averaged AM results are 
tabulated in the EGA columns of Table |^ The results are given as a fraction of the vr 
pair current contribution, as this is the dominant term. These results are compared to full 
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TABLE IV. A comparison of anapole moment estimates from a one-body Fermi gas average 
with full two-body shell model results. DDH best-value couplings are used, and no short-range 
correlation function corrections are included in either set of results. The labels PC and PNC denote 
whether the nucleon absorbing the photon has a PC or PNC meson- nucleon coupling. 



ID/ex 10^ 


vr pair 


p pair PC 


p pair PNC 


Lj pair PC 


uj pair PNC 


7r7r7 


PPl 


p-n-jpc 


FGA 


110 


13.0% 


-19.0% 


-0.4% 


8.1% 


-34.9% 


6.6% 


0.5% 


SM 


67 


12.9% 


-18.2% 




8.6% 


-24.0% 






FGA 


-75 


12.8% 


-18.2% 


-0.3% 


7.8% 


-35.5% 


7.8% 


0.0% 


SM 


-27 


15.4% 


-21.5% 




12.8% 


-29.4% 







two-body SM results, similarly normalized to the SM vr pair current AM value. The absolute 
vr pair current results are also given for both calculations. 

We see from the table that, while the Fermi gas average tends to overestimate the AM 
contribution by a factor of ~ 2-3, compared to the SM, the Fermi gas and SM agree very 
well on the relative values of the various contributions. (The comparison is less impressive 
for Tl than for Cs, but the Fermi gas parameters used for both nuclei were tailored to Cs.) 
This suggests that the one-body average AM values should be reliable indicators of which 
exchange current contributions are important. 

The Fermi gas model is an independent particle model. The SM, while incorporating 
certain correlations, omits the high-momentum components of the Hilbert space necessary 
for describing the short-range hard core. While the SM (and associated Fermi gas) short- 
comings could in principle be corrected by introducing effective operators and wave function 
renormalizations, in practice this is never done. Instead, most frequently the omitted short- 
range physics is mocked up by a correlation function which, in SM PNC studies, is often 
taken from Miller and Spencer [Q, 



/(n2) = l-(l-ferf2)e-'^'--, (61) 

with a = 1.1 /m~^ and b = 0.68 /m~^. This correlation function reduces two-body matrix 
elements by ~ 25-30% for vr currents, 75-80% for p and uj currents, ~ 80% for rni currents, 
and ~ 90-95% for pp and pn. 

No short-range correlation corrections have been included in the results of Table 
It is thus apparent that the true j'^J'pc , (the most complicated current), and jp^'^ 
exchange current contributions (with short-range correlations included) would be ~ 1% of 
the dominant vr pair result. It is then reasonable to ignore these unimportant but complicated 
exchange currents, evaluating all others with the full two-body SM density matrix, modified 
by the Miller-Spencer correlation function. While a complete list of the two-body AM 
operators is too long to list here, the dominant vr operator is found to be 

a^-^ = '";;^;'^ (r(l)-r(2)-.3(l)r3(2)) 



X 



xla{l) + xld{2) + V2^{xl[Y2{n^) a(l)]i + a;2[F2(fi2) ® ^(2)]i} 
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--(a(l)-Vi-a(2)-V2) 



X ((x^ + xl)x + V2^{xl[Y2{ni) f]i + xl[Y2{n2) ® f]i} + |— f 



2 m„ 



a; 



(62) 



where x = xi — X2- The numerical results for the sum of all exchange current contributions 
to the Cs and Tl AMs is given in Table [VT] . 



C. Nuclear Polarization Contributions 

As illustrated in Fig. ^ the two-body PNC A^A^ potential perturbs the ground state, 
mixing it with excited states of opposite parity. The resulting odd-parity ground state 
component allows the ordinary (vector) El current to couple to the ground state. The 
first-order perturbation theory AM is thus 

^ {I\a^\n){n\H^;lc\I) , {I\H^plM{n\a''\I) 

where |/) is the unperturbed ground state of good parity and the sum extends over a 
complete set of nuclear states n of angular momentum I and opposite parity. The operator 
aX is obtained by plugging the ordinary electromagnetic current into Eq. |2D|, 



+ \ ® [/is + /i^rs {i)] | (64) 

where Hs = 0.88 and = 4.706. 

The summation over a complete set of intermediate SM states for ^^^Cs or 205^| jg 
impractical either directly or by the summation-of-moments method discussed in Ref. 
and below. However, because no nonzero El transition exists among the valence orbits 
(e.g., the /iii/2 and (77/2 orbitals have opposite parity but cannot be connected by a dipole 
operator), an alternative of completing the sum by closure, after replacing l/AE^ by an 
average value (l/AE) is quite attractive 

_ ^ {I\a''\n){n\H^;l c\I) + (n|a^|/) 



AE„ 



The resulting product of and Hpj^Q contracts to a two-body operator, so that only the 
two-body ground state density matrix is needed, a considerable simplification. (No three- 
body terms arise because the absence of El valence transitions guarantees they vanish for 
our SM spaces.) 

The closure approximation can be considered as an identity, clearly, if one knows the 
correct (l/E), that is, how to parameterize the relationship between the weighted and 
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non-energy-weighted sums. In practical terms, this means demonstrating that a systematic 
relationship exists between (l/E) and some experimentally known quantity, such as the 
position of the El giant resonance. Note that the El operator is closely related to the 
anapole operator . 

To investigate the systematics we completed a series of exact calculations in Ip— and light 
2sl(i— shell nuclei (''Li, ^^B, 17,19, 2ip^ 2i,23]v^g^'j^ evaluating both the (l/E) and non-energy- 
weighted sums. First, the ground states are determined from full Ohu diagonalizations. The 
polarization sum involves the complete set of Ihu states that connect to the ground state 
through the anapole operator. The summation was performed by exploiting a variation of 
the Lanczos algorithm to evaluate the effect of the nuclear propagator 1/Eg g — H (see Sec. 
V.D). The algorithm efficiently completes the sum via moments, even though the dimensions 
of the Ihuj bases ranged up to ~ 500,000. The appropriate closure energies were found not 
only for the anapole polarization sum, but also for El operator. This allowed us to compare 
the (l/E) appropriate for the AM calculation with that appropriate for photoexcitation. As 
photoexcitation response functions have been mapped in many nuclei, this in turn allows us 
to relate the anapole {l/E) to an experimental observable. 

The results show that the anapole and photoexcitation average excitation energies track 
each other very well, provided one takes into account the three isospins contributing to 
Hpl;(j. Measured as a fraction of the l/i?-weighted giant dipole average excitation energy, 
which is (l/E)^^ ~ (22-26) MeV for these nuclei, the appropriate effective energies for the 
anapole closure approximation are 0.604 ± 0.056 for and (isoscalar channel), 0.899 ± 
0.090 for (isovector channel), and 1.28±0.14 for /i^ (isotensor channel). The larger (l/E) 
for h^p and enhances the isoscalar contribution to the anapole polarizability. The small 
variation in (l/E), once the isospin dependence is recognized, supports the notion that we 
can connected the closure result to the true polarization sum. 

Inspired by the nuclear systematics we found above, we estimate T=0,l,2 closure energies 
from known El distributions, that is, we fix the anapole closure energy as 0.6, 0.9, and 1.28 
of the El closure energy evaluated from the experimental dipole distribution. For ^^^Cs 
41|, this gives 9.5, 14.1, and 20.2 MeV, respectively. The corresponding ^os^jj values are 
8.7, 12.9, and 18.5 MeV. The ground-state expectation values for the contracted two-body 
effective operator {a^, Hplj^^} are then evaluated from the SM two-body density matrices 
for Cs and Tl. The Miller-Spencer correlation function is again included in the two-nucleon 
matrix elements of Hpj^^. The resulting polarization contributions are given in Table VI. 



V. EXPERIMENTAL CONSTRAINTS, RESULTS, AND UNCERTAINTIES 

In this sections we discuss atomic PNC experiments that determined (or limited) the 
AMs of ^^'^Cs and ^°^T1, other experimental tests of the PNC hadronic weak interaction, 
and the consistency of the AM results with these other tests. We also discuss nuclear 
structure uncertainties in the interpretation of the AM measurements. 
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TABLE V. Nuclear systcmatics found in light odd-proton nuclei: The second column shows the 
functional dependences of SM results for the direct anapole polarization sums and the third column 
shows the forms for sums by closure approximation using the closure energy {1/E)ei which is 
derived from the weighted El sum rule (also evaluated in the SM). The same normalization has 
been applied to the second and third columns. By comparing these two columns, it is apparent that 
in order for the closure approximation to be correct, the anapole closure energies (l/-E-)yijv/(ji=o,i,2) 
should be different from {1/E)ei- In columns 4-6 we express (l/^)^7v/(ji=o i 2) units of 
thus a value less than one means that the appropriate anapole average excitation energy is lower 
than the corresponding average over the photoexcitation peak. Note the closure result faithfully 
reproduces the correct — h^, combination. We omit the dependence on and h^, because the net 
isovector contribution is almost entirely from Z^. In the case of ^^F, the lowest, nearly degenerate 
l/2~ state was removed from all sums. 



Nucleus Direct Pol. Sum Closure with {1/E)ei (V-^)am(o) (V-^)am(i) (^/^^aI 



U - 0.34(/i" -I- 0.58h!i) + 0.05W 0.80/^ - 0.20(W -|- OSSh)^) + 0.05/i^ 0.59 0.80 To 

"B /,r - 0.53(/i° -l-0.52/iO) -|-0.05ft^ 0.89/,r - 0.37{/i° -|- 0.52/i° ) -l- 0.07/1^ 0.70 0.89 1.4 

-0.60{h°p + 0A8hl) + 0.0Ah^ 1.02/„ - 0.40(ft0 + 0.46?i° ) + 0.05/i^ 0.66 1.02 1.2 

i^F /,r - 0.33(/i° -1-0.56/1°) + 0.02/i2 0.90/^ - 0.19(/i° + 0.59/i° ) + 0.03/1^ 0.58 0.90 1.5 

^^F 0.41(/i0 -1-0.55/1°) + 0.03/1^ 0.97/^ - 0.24(/i° + 0.54/1° ) + 0.04/i2 o.60 0.97 1.3 

2iNa - 0.57(/i° -f 0.51/1°) + 0.02/i2 0.77/^ - 0.31(/i° + 0.49/i° ) + 0.03/i? 0.54 0.77 1.5 

23Na - 0.67(/i° -1-0.53/1°) -1-0.05/1^ 0.95/^ - 0.38(/i° -|- 0.52/i° ) -|- 0.07/i| 0.57 0.95 1.4 



TABLE VI. Decomposition of the SM estimates of the anapole matrix element (/||^i||/)/e 
into its weak coupling contributions. 



Nucleus 


Source 


fn 




hi 


hi 




hi 


133Cs 


nucleonic 


0.59 


0.87 


0.90 


0.36 


0.28 


0.29 




ex. cur. 


8.58 


0.02 


0.11 


0.06 


-0.57 


-0.57 




polariz. 


51.57 


-16.67 


-4.88 


-0.06 


-9.79 


-4.59 




total 


60.74 


-15.78 


-3.87 


0.36 


-10.09 


-4.87 


205 rpj 


nucleonic 


-0.63 


-0.86 


-0.96 


-0.35 


-0.29 


-0.29 




ex. cur. 


-3.54 


-0.01 


-0.06 


-0.03 


0.28 


0.28 




polariz. 


-13.86 


4.63 


1.34 


0.08 


2.77 


1.27 




total 


-18.03 


3.76 


0.33 


-0.30 


2.76 


1.26 
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FIG. 7. Atomic parity mixing induced by exchange. 

A. Constraints from the Nuclear Anapole Moments of ^'^^Cs and ^os-pj 



A thirty-year program to study atomic PNC has yielded in the past few years 
exquisitely precise (sub 1%) results. The primary focus of these studies has been to obtain 
accurate values of the strength of direct exchange between electrons and the nucleus. The 
PNC effects are dominated by the exchange involving an axial Z° coupling to the electron 
and a vector coupling to the nucleus. The nuclear coupling is thus coherent, proportional to 
the weak vector charge, Qw ~ Z{\ — 4 sin^ 6*^/) — ~ — A^, and independent of the nuclear 
spin direction. It is widely recognized that these atomic measurements are important tests of 
the standard electroweak model and its possible extensions, complementing what has been 
learned at high energy accelerators that directly probe physics near the pole [^,0 ■ 



In heavy atoms the weak electron-nucleus interaction will induce a small P— wave parity 
admixture in an atomic S orbital on the order of parts in 10^^. This will produce, in a 
transition that is normally Ml, a small E\ component. The PNC signal will be easier to 
detect if the parity-allowed Ml transition is hindered, as the observable depends on the 
E\IM\ ratio. The forbidden Ml transitions of 6S'i/2 — >■ 75*1/2 in Cs and 6P1/2 — ^ 7Pi/2 in Tl 
are two examples of this sort. Moreover, the structure of these atoms is comparatively simple, 
allowing theorists to extract the underlying weak couplings from the PNC observables. 

One popular atomic technique exploits the linear Stark response to an applied static elec- 
tric field. A coordinate system in the atom is established by mutually perpendicular Stark, 
magnetic (for producing the Zeeman spectrum of states that can be populated by optical 
pumping), and laser (stimulating the El transition) fields. The "parity transformation" is 
accomplished by inverting these fields. The PNC signal is associated with any difference 
seen in the interference between the Stark, PNC i?l, and hindered Ml amplitudes after 
various reversals of the coordinate system. The elimination of spurious signals associated 
with imperfect field reversals and other sources of systematic error is a tedious task. A 



recent review of the Cs and Tl experiments can be found in [45 



The dominant axial(electron)-vector(nucleus) atomic PNC interaction is independent 
of the nuclear spin (see Fig. There is also a tree- level contribution to atomic PNC 
that is nuclear-spin-dependent, where the Z° exchange is vector(electron)-axial(nucleus). 
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This contribution is highly suppressed because the vector electron weak coupling is small, 
Qy^ = —(1 — 4sin^6'vK) ~ —0.1, and the nuclear coupling is no longer coherent. But, given 
sufficiently accurate (< 1%) measurements, this suppressed signal can be cleanly extracted 
by studying the hyperfine (and thus nuclear spin) dependence of the PNC measurements. 

In Sec. II we noted that the nuclear AM will also generate a nuclear-spin-dependent 
weak interaction between the electron and the nucleus, thus contributing in combination 
with tree-level V(electron)-A(nucleus) exchange. Furthermore other 0{GF<y) radiative 
corrections also contribute to that spin dependence, with the hyperfine interaction between 
the electron and nucleus (see Fig. ||) of particular importance because of the coherent 
coupling. While the naive expectation is that radiative corrections will indeed be corrections 
of strenth ~ a relative to the tree-level contribution, the small vector coupling of the Z° to 
the electron combined with the A^^^ growth of the anapole moment leads to a surprise. The 
AM becomes the dominant source of nuclear-spin-dependent atomic PNC for A > 20 . 
This guarantees not only that the nuclear spin dependence is signifcant for heavy atoms, 
but also that the AM contribution might be deduced from the measurements. 

The nuclear-spin-dependent (NSD) PNC electron-nucleus contact interaction which gen- 
erates the parity mixing can be expressed as 

HpNE = ^^tota-Ipir), (66) 

f^tot = tz" + f^hf + (67) 

— * 

where / and p(r) are the nuclear spin and density, a the usual Dirac matrix of the electron, 
and K a dimensionless constant which characterizes the strength of the PNC. (Note that 
our definition of k is different from the one given by Khriplovich and others by a factor 
(— 1)'^"*"^/^"*"'(/ + l/2)/(/(J + 1)), where / is a single-particle orbital angular momentum. The 
Khriplovich definition thus assumes a single-particle picture, though there are examples of 
nuclei where the dominant single-particle orbital is characterized by an / that is naively 
inconsistent with the many-body /, e.g., I ^ I ± 1/2.) The k, subscripts denote contribu- 
tions from Z^ exchange, the hyperfine interaction correction, and the AM. From the ^^^Cs 
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(extracted by Flambaum and Murray |T9|) and ^05^] results one finds 

KtotC^^Cs) = 0.112 ±0.016, 
KtotC^^Tl) = 0.29 ± 0.40 Seattle, 

/ttot(2°^Tl) = -0.08 ±0.40 Oxford. (68) 

Henceforth we will focus on the Seattle Tl result, as this proves to be more restrictive 
than the Oxford result in the parameter space of PNC hadronic couplings favored by other 
experiments. (The Oxford AM result is quoted with opposite signs in different sections of 

and the accuracy of the spin-independent measurement is considerably less than that 
of the corresponding Seattle measurement. These observations contributed to our decision 
to focus on the result of 0.) We treat the Tl constraint as one on the principal isotope 
205rp| (70.5%). The other stable isotope, ^os^j] (29.5%), differs in structure only by a pair of 
neutrons, and thus should have very similar properties. 

The contribution is 

, . 2n ^ allEf=l^(^)r3(^)||/) 
Kzo = -y(1 -4sm Ow) ' ^ 

with the axial vector coupling = 1.267 and sin^ 9\y = 0.2230. Here 1 1 denotes a matrix ele- 
ment reduced in angular momentum. The reduced matrix element of / is ^J(J ± 1)(2J ± 1). 
The Gamow- Teller matrix elements, taken from the SM studies, are —2.305 (^^^Cs) and 
2.282 (^°^T1), not too different from the corresponding single-particle (s.p.) values of —2.494 
(unpaired 1(77/2 proton) and 2.449 (3si/2 proton). This yields: 

KzoC^^Cs) = O.OUO, (70) 
k^o(205t1) = -0.127. (71) 

Note that the inclusion of one-loop standard model electroweak radiative corrections modify 
these results, reducing the isovector contribution substantially and inducing a small isoscalar 
component. 

For the hyperfine correction, from the measured nuclear weak charge and magnetic mo- 
ment, Bouchiat and Piketty [|11|] find 

KhfC^^Cs) = 0.0078, (72) 
KhfC^^Tl) = 0.044. (73) 



Note that the conversion of the notation of Ref. iTTH to ours is 



' ' mm 

By subtracting k^o and Hhf from Utot we obtain the AM contribution 



(74) 



kam(^^^Cs) = 0.090 ± 0.016, (75) 
kam(^°^T1) = 0.376 ± 0.400. (76) 
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These values are related to the nuclear AMs by 



47rav/2(/||a||/)/e 
'^^^^ - G,M- (/||/||J) ' 

where a is the anapole operator. As our results for (/||a||/)/e are expressed in terms of the 



PNC meson-nucleon couplings in Table VII, we have the needed AM coupling constraints. 



B. Constraints from Nuclear PNC Experiments 

The nuclear experiments measuring an interference between PC and PNC amplitudes 
generally fall into four types 

1. Measurement of the longitudinal asymmetry in a scattering experiment (e.g., pp, 
pd, or pa) . 

2. Measurement of the circular polarization of photons emitted in a nuclear decay 
(e.g., ^*F, ^-^Ne) or reaction (e.g., np — > ^7). 

3. Measurement of the asymmetry of photons emitted in the decay of a polarized 
nucleus (e.g., ^^F) or in a polarized nuclear reaction (e.g., np d'j, nd — > ^7). 

4. Measurement of the degree of spin rotation for polarized neutrons through various 
targets (e.g., p, (i,''He). 

It is unfortunate that only a single A^A^ PNC scattering observable, the longitudinal ana- 
lyzing power Al for p+p, has been successful p6|-^. (Experiments have been done at 13.6, 
45, and 221 MeV.) These results have been supplemented by a number of PNC measure- 
ments in nuclear systems, where accidental degeneracies between pairs of opposite-parity 
states can produce, in some cases, large enhancements in the PNC signal. Unfortunately 
not all of these results are readily interprettable because of nuclear structure uncertainties. 



Those that can be analyzed with confidence p[ include A^ for p + a at 46 MeV p9[, the 



circular polarization of the 7-ray emitted from the 1081 keV state in ^^F |5^, and A^ for 



the decay of the 110 keV state in polarized ^^F [ET|. These examples involve either few-body 



systems, where quasi-exact structure calculations can be done, or special nuclei in which the 
PNC mixing matrix elements can be calibrated from axial-charge (3 decay [^. An analysis 
of these results, which have been in hand for some time, suggests that the isoscalar PNC 
A^A^ interaction - which is dominated by p and uj exchange - is comparable to or slightly 
stronger than the DDH "best value," whereas the isovector interaction - dominated by vr 
exchange - is significantly weaker (< 1/3) [^. Because one expects the isovector channel to 
be governed by neutral currents and to receive potentially significant light sea-quark con- 
tributions, there is considerable interest in testing this result. The Cs and Tl AM results 
provide one possible cross check. 
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TABLE VII. PNC observables and corresponding theoretical predictions, decomposed into the 
designated weak-couphng combinations. 



Observable 


Exp.(xlO'^) 


/^-0.12/ii -0.18/ii 


h% + 0.7/1° 


h\ 




"■w 


hi 


J TIT) 1 ^ ^ \ 

^^'^(13. 6) 


-0.93 ± 0.21 




0.043 


0.043 


0.017 


0.009 


0.039 


(45) 


-1.57 ± 0.23 




0.079 


0.079 


0.032 


0.018 


0.073 


(221) 


0.84 ± 0.34 




-0.030 


-0.030 


-0.012 


0.021 




^r(46) 


-3.34 ± 0.93 


-0.340 


0.140 


0.006 




-0.039 


-0.002 


P,(18f) 


1200 ± 3860 


4385 
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-44 


A,(19f) 


-740 ± 190 


-94.2 


34.1 


-1.1 




-4.5 


-0.1 


(||Ai||)/e, Cs 


800 ± 140 


60.7 


-15.8 


3.4 


0.4 


1.0 


6.1 


(ll^i||)/e, Tl 


370 ± 390 


-18.0 


3.8 


-1.8 


-0.3 


0.1 


-2.0 



C. Results 

The constraints on PNC meson-nucleon couplings of Table |V11| are displayed graphically 



in Fig. 1^. Although there are six independent couplings, two combinations of these, one 
isoscalar and one isovector, dominate the observables: /^r — 0.12/i^ — 0.18/i(^ and -|- 0.7/i° . 
The decomposition of Table |V11] thus uses these two degrees of freedom along with h?^ and 
the residual contributions in /i^, and h}^. The \o error bands of Fig. |^ are generated from 
the experimental uncertainties, broadened somewhat by allowing uncorrelated variations in 
each of the four minor degrees of freedom (that is, and the residuals in in /i^, , and /i^) 
over the DDH broad "reasonable ranges." Note that only a fraction of the region allowed by 
the Seattle Tl constraint is shown: the total width of the Tl band is an order of magnitude 
broader than the width of the Cs allowed band, with most of the Tl allowed region lying 
outside the DDH "reasonable ranges" (i.e., in the region of negative /vr — 0.12/;,^ — O.lS/i^ 
and positive /i^ + 0.7/i°). That is, the bulk of the Seattle Tl band corresponds to an AM 
value opposite in sign to that expected theoretically, given what we know experimentally 
about PNC meson-nucleon couplings. The corresponding Oxford Tl band (not illustrated) 
includes almost all of the parameter space in Fig. |^, as well as a substantial region outside 
the bounds of the figure, to the lower left. 

The weak coupling ranges covered by Fig. ^ correspond roughly to the DDH broad 
"reasonable ranges." Thus the anapole constraints are not inconsistent with the theoretical 
"ball-park" estimates. However, the detailed lack of consistency among the various mea- 
surements is disconcerting. Before the anapole results are included, the indicated solution 
is a small /tt and an isoscalar coupling somewhat larger than, but consistent with, the DDH 
best value, — (/ip + 0.7/;.°)^^^ ~ 12.7. But the AM results agree poorly with this solution, 
as well as with each other. In particular, the precise result for ^^^Cs tests a combination of 
PNC couplings quite similar to those measured in A^i^'^Y') and in A^", but requires larger 
values for the weak couplings. 

Despite substantial differences between our work and that of Flambaum and Murray 
[T9| , the predicted AMs from these two calculations are in relatively good agreement. The 



corresponding interpretations, however, are quite different. Flambaum and Murray adopted 
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the viewpoint that the Cs AM result could be accommodated by a value /,r ~ 9.5, about 
twice the DDH best value, /^f^ ~ 4.6. (The DDH reasonable range is 0-11.4, in units of 
10~^.) The difficulty with this suggestion is its inconsistency with P^(^^F), a measurement 
that has been performed by five groups. The constraint from this measurement is almost 
devoid of theoretical uncertainty 

- 0.6 <U- 0.11/iJ - 0.19/ii < 1.2. (78) 

If one allows and hlj to vary throughout their DDH reasonable ranges, one finds — 1.0 < 
fjr < 1.1, clearly ruling out ~ 9. There is also some tension between the Cs band and 
those for p + a and A^{^^F). 

Thus, unfortunately, the hint of a consistent pattern of weak meson-nucleon couplings 
that was emerging from nuclear tests of the weak hadronic current is disturbed when the Cs 
and Tl results are added. 



D. Operator Renormalization and other Nuclear Structure Issues 

It thus appears that the calculated value of the Cs AM, using weak meson-nucleon 
couplings determined from A^A^ and nuclear experiments, is significantly smaller than the 
measured value. While there are several questions that could be raised about this conclusion, 
perhaps the most difficult one is the quality of the nuclear structure calculations for Cs 
and Tl: what error bar should we assign because of the inherent uncertainties in such 
calculations? 

Despite the rather extensive theoretical literature on AMs, it would be fair to characterize 
the general quality of the associated nuclear structure work as unsophisticated. Much of the 
previous work is based on extreme single particle models and employs effective one-body 
PNC potentials, a choice that tends to obscure the discrepancies apparent in Fig. ^. Only a 
few attempts have been made to estimate the effects of correlations, even in schematic ways. 



In [|Tl| quenching factors were introduced as a phenomenological correction to single-particle 



estimates. Solid motivation for this approach can be found in classic studies of magnetic 



moments and Gamow- Teller transitions in nuclear physics. In |]15[ single-particle calculations 
were corrected for core polarization effects, employing a realistic g-matrix interaction but a 
very simple set of particle-hole excitations. Despite the highly truncated model space, this 



may be the only paper, other than our work here and in earlier papers P,pi][ , to use a realistic 



interaction in calculations of the Cs and Tl AMs. Finally, in Ref. |T3| core polarization 
effects were evaluated in the random phase approximation, but with a schematic zero-range 
spin-spin residual interaction. 

One factor limiting what can be done is the challenge of completing the polarization 
sum: apart from the work referenced above performed this sum state- by-state. Such 

a summation technique rules out a sophisticated ground-state wave function: the number 
of opposite-parity eigenstates connecting to the ground state by the El operator would 
be enormous. The two attempts to move beyond direct summation have come from our 
studies. In summation to a complete set of Ihu states for ^^F was carried out by a 
Lanczos algorithm moments method. In this approach one recognizes that the quantity of 
interest is the distribution of the vector a^\I) over the full set of Ihu eigenstates: if that 
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FIG. 9. Constraints on the PNC meson couplings (xlO ) that follow from the results in Table 
VII. The error bands are one standard deviation. The pp band is the union of 13.6, 45, and 



221MeV results. 
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distribution is known, it can be weighted by l/{Eg,s. — En) and dotted with {I\Hpy to 
generate the polarization sum. Instead of diagonahzing a very large matrix of dimension N , 
where N is the number ITiuj eigenstates, to get the eigenvalues En and eigenstates needed 
to do this sum state-by-state, the Lanczos method maps the large matrix into a series of 
smaller matrices of dimension A^' = 1,2,3..., where A^' -C A^. This mapping extracts exact 
information from the original large matrix, the 2N' — 1 lowest moments of the vector a^| J) 
over the Ihuo eigenspectrum. It is readily seen that the distribution must be very well 
determined after a modest number of iterations, A^' ~ 50. There is a variation of this 
algorithm that uses the information in the Lanczos matrix to construct the effect of the 
Green's function [^]: it is obvious physically that one can obtain the Green's function from 
the detailed moments construction. (The algorithm develops the Green's function acting 
on a vector as an expansion in the Lanczos vectors, with the the coefficients of the vectors 
updated with each iteration |^3|. The method is thus exact in a numerical sense, allowing 
one to evaluate the convergence.) This was the method used in the present study of p— 
and sci— shell nuclei, to assess average excitation energies. We have applied this method 
in cases where A^ ~ 10^, and it is possible with modern machines to tackle problems of 
dimension ~ 10^ in this way. Unfortunately, given the complexity of our ^^^Cs ground state 
wave function, the dimension of the negative-parity space required to saturate the El sum 
is substantially larger than 10®. Thus this technique, while exceedingly powerful, cannot be 
applied to a case like ^^^Cs, at least at the present time. 

Because we felt it was important to use a realistic large-scale SM wave function in 
describing the ^^^Cs ground state, another method was needed to evaluate the polarization 
sum. We did this by closure, which was tractable in part because of an attractive property of 
the canonical ^^'^Cs SM space, no nonzero matrix elements of oX . In our view there are two 
worrisome features of this calculation. The first is the reliability of the average excitation 
energy estimate, which we defined as the ratio of the non-energy-weighted to l/i?„-weighted 
sums. We performed a large set of calculations in lighter nuclei, using the exact Lanczos 
Green's function method described above, to calibrate the method. The average excitation 
energies, normalized to the photoexcitation El peak and evaluated for each isospin channel, 
proved to be very stable. One cannot prove that the extrapolation to heavy nuclei like Cs 
and Tl is valid, clearly: perhaps there is some systematic evolution with neutron excess. On 
the other hand, the naive expectation is that the method should improve with A, as the El 
profile tends to become more collective in heavier nuclei, and as the spin-orbit force tends to 
remove El strength from low excitations: the closure approximation is clearly exact in the 
limit of an infinitely narrow El resonance. Because the measured Cs AM is large, one would 
need a substantial amount of strength quite low in the Cs spectrum to enhance the 1/ En sum 
and thus "fix" the SM calculation: this is unexpected and, while the aX and photoexcitation 
El operators are somewhat different, there is no evidence in the photoexcitation distribution 
for such strength . 



The second question is the adequacy of our ground-state wave function: though the Cs 
and Tl SM calculations are serious efforts, numerical limitations forced restrictions on the 
proton and neutron occupation numbers. The unrestricted lg-j/2 — '2d^/2 — 3si/2 — 2(^3/2 — 
l/iii/2 SM calculation was not attempted. Furthermore, it is well known that even full-shell 
calculations often must be renormalized phenomenologically. Two operators closely related 
to the AM, the Gamow- Teller and Ml operators, are well-studied examples [Bl]. In Table 
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TABLE VIII. Magnetic moments of ^^^Cs and ^"^Tl measured in nuclear magnetons. 





s.p. 


SM 


exp. 


133Cs 


1.72 


1.65 


2.58 


205 -p^ 


2.79 


2.58 


1.64 



Villi our Cs and Tl SM magnetic moment values are compared to the experimental and 
s.p. values. The SM and uncorrelated s.p. values are not that different, and both differ 
significantly from experiment. The conclusion is that potential important physics is absent 
in our truncated SM calculations. 

The deviations of magnetic moments from the Schmidt line (or s.p. values) around the 



Pb region have been extensively studied by Arima et al. |54]. The deviations from the s.p. 



predictions can be described as a set of corrections to the bare gyromagnetic factors 

{iWmil^N ={\ + SgfMms.p. + {\ + 59\'^)mr,\\)s.p. + (O.SS + 6gf)m\)s.r>. 

+ (4.70 + 5<7«)(||ar3||)..,., (79) 

These factors represent the operator and wave function normalization corrections that would 
result from a faithful treatment of the omitted parts of the Hilbert space. Equivalently (and 
perhaps more appropriately) one can quote this result in terms of renormalized matrix 
elements 

(J||/i||/)//i;v = ^(/||/l|/)ren + ^(/||rr3||/).e„ + 0.88(/||a||/),en + 4.70(/||ar3||/),^^ (80) 



The fit of [54] gives the following quenching for the spin matrix elements near Pb 



(/||a||/),en = 0.86(1 (81) 
(/||(Tr3||/),en = 0.54(1 larsll),,^,. (82) 

Although there exists no such large body of data on the anapole moment operator, 
we now explore whether some tentative conclusions can be drawn about effects of missing 
correlations on that operator. We begin with the observation that the effects of correlations 
on a many-body operator are expected to be quite similar to their effects on the one-body 
equivalent of that operator. (One specific illustration of this is detailed in |§2[-) Thus we 
start by looking for the one-body equivalent of the anapole polarization operator. The most 
general spin-isospin form for a rank-one operator is 

= ^(«f°^r+apV3r+afa + a«r3a + af [r2®a]i + 4'V3[r2®a]i). (83) 

As the average excitation energy is measured in units of hoo, the bare couplings are 
dimensionless. We then evaluate matrix elements of this one-body operator and of the full 
polarization sum (chosing DDH "best-value" meson-nucleon couplings) in a single-particle 
model for a variety of nuclei in the Pb and Sn regions, fitting the coefficients of the one-body 
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operator to reproduce the polarization results. The results for Tl (Pb region) are presented 



in a series of three tables, Tables |X|, 0, and ^I|, giving, respectively, the comparison of 



the calculated s.p. polarization results with those generated by the effective operator, the 
best fit values found for the coefficients of the effective operator, and the matrix elements 
of the various terms in the effective one-body operator. The following three tables give the 
analogous results for Cs (Sn region). Calculations were done with no spin-orbit potential as 
well as with a spin-orbit potential of strength -O.lhuja ■ I: the results show little sensitivity 
to the spin-orbit contribution. 

The tables show that the orbital contributions to the effective operator are neglible: the 
dominant terms are the spin and spin-tensor operators, with the former (folding the results 
of Tables ^ and X] and of Tables X.111 and XIV ) accounting typically for about 70% of the 
AM strength. Furthermore the spin isoscalar and spin isovector operators contribute with 
the same relative sign, with the isovector contribution larger. It follows for ^°^T1, where the 
single particle assignment is 3si/2, eliminating both the spin-tensor and orbital contributions, 
that the effective AM operator is very similar to the magnetic moment operator, and thus 
should be renormalized in a very similar way. From Table Vlll one concludes that our SM 
estimates are not sufficiently quenched, overestimating the Tl AM by about a factor 1.6. 
The consequence of this would be to broaden the allowed Tl band (only partially shown) in 
Fig. 1^ proportionately. 

The case of ^^^Cs is more difficult in that the spin-tensor operator now plays a significant 
role: the s.p. assignment is 1(77/2- This operator does not arise as a bare operator in 
Gamow- Teller, Ml, or other familiar responses. Our approach is somewhat unsatisfactory, 
but perhaps of some help. In Table pCV| we compare s.p. and full Ip— and 2sl(i— shell SM 
calculations of magnetic moments with the experimental values for a series of light nuclei. 
This seems to establish that, in these nuclei, the bulk of the needed renormalization of s.p. 
estimates does come from the SM (sweeping under the rug issues like exchange currents, etc.). 
In Table pCVl| we make a similar comparison of s.p. and SM AM operator matrix elements. 
The pattern of significant quenching of spin matrix elements again emerges from this purely 
theoretical comparison. In the case of the spin-tensor operator, the renormalizations do not 
seem very large, nor do they appear to follow a simple pattern. While there are cases of 
modest spin-tensor matrix element enhancement when the full-shell correlations are turned 
on, these enhancements are smaller than the quenching that occurs in the spin matrix 
elements. The overall tendancy of the correlations is to suppress the AM prediction. 

While these arguments are of a hand-waving nature, they favor the conclusion that bet- 
ter SM calculations will produce a somewhat smaller, not larger, predicted Cs AM. The 
dominant missing physics appears to be insufficient quenching of the spin matrix elements. 
This will clearly exacerbate the discrepancies apparent in Fig. As a full-shell calcula- 
tion for ^^^Cs will likely become feasible within the next few years, there may soon be an 
opportunity to demonstrate that improved calculations will produce a smaller AM. 



VI. CONCLUSIONS 

Recent atomic PNC measurements in ^^'^Cs reached a new level of precision that led, for 
the first time, to detection of the hyperfine dependence of the signal. New measurements in 
Tl have also imposed important constraints on nuclear-spin-dependent atomic PNC. This 
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TABLE IX. Comparison of calculated s.p. polarization anapole moments (||a||)/e in the ^'^^Pb 
region with results for a fitted phenomenological effective operator. 



Nucleus 




No s.o. 




With s.o. 






Calc. 




Fit 


Calc. 


Fit 



207t1(3s7') 
207t1(24^) 
207Tl(2(i^) 
207pb(3p7^) 

207pb(2/J-i) 

20961(2//) 
209Bi(2/J) 



209 T 



^Bi(l/l9) 

209pb(2/9) 



-578 

759 
-780 
-131 

161 
-180 

970 
-919 
1154 

224 



-593 

763 
-889 
-122 

154 
-190 

924 
-1012 
1198 

232 



-542 
699 

-691 

-132 
158 

-184 
881 

-821 
990 
212 



-536 

692 
-825 
-123 

151 
-194 

830 
-949 
1057 

220 



TABLE X. The fitted parameters in ^ospb region. 



dj^ Qjj^ dg Qig Ctp dp 

no s.o. 0.990 1.458 -95.838 -146.159 -243.094 -366.696 

with s.o. -0.721 0.432 -84.580 -134.308 -224.986 -348.570 



TABLE XI. The single particle reduced matrix elements used for ^^Spb region fits. 

I T3I ^ T3^ [I2 ® T3[Y2 (8) d]l 

20'^Tl(3s7^) 0.000 0.000 2.449 2.449 0.000 0.000 

20'^Tl(2d3^) 4.648 4.648 -1.549 -1.549 -0.618 -0.618 

20'^Tl(2dB^) 5.797 5.797 2.898 2.898 0.330 0.330 

20'^Pb(3p7^) 1.633 -1.633 -0.816 0.816 -0.651 0.651 

20'^Pb(3p3^) 2.582 -2.582 2.582 -2.582 0.206 -0.206 

20'^Pb(2/7^) 8.281 -8.281 -2.070 2.070 -0.661 0.661 

209Bi(2/5') 8.281 8.281 -2.070 -2.070 -0.661 -0.661 

209Bi(2/7) 9.621 9.621 3.207 3.207 0.426 0.426 

209Bi(l/i9) 17.162 17.162 -2.860 -2.860 -0.761 -0.761 

209pb(2/9) 13.984 -13.984 3.496 -3.496 0.507 -0.507 
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TABLE XII. Comparison of calculated s.p. polarization anapole moments 
region with results for a fitted phenomenological effective operator. 



a\\)le in the ^^SSn 



Nucleus 




No s.o. 




With s.o. 






Calc. 




Fit 


Calc. 


Fit 



I3iln(2p7') 

I3iln(2p^) 

i3iSn(3/7') 
i3iSn(2d^) 
i3iSn(24i) 

i33Sb(l5j) 
i33sb(2(il) 
i33Sn(2/J) 
i33Sn(l/i9) 



433 
641 

-484 
103 

-124 
137 
788 

-610 
169 

-169 



422 
644 

-450 
93 

-109 
144 
751 

-534 
168 

-159 



409 
567 

-440 
102 

-125 
127 
684 

-549 
158 

-171 



429 
636 

-449 
89 

-106 
136 
736 

-539 
157 

-160 



TABLE XIII. The fitted parameters in ^^Sgn region. 



tt^ Qij^ Qjg Qig dp dp 

no s.o. 3.284 2.327 -52.990 -0.790 -182.832 -271.131 

with s.o. 1.807 1.315 -44.608 -80.981 -176.206 -260.346 



TABLE XIV. The single particle reduced matrix elements used for ^^^Sn region fits. 

I T3I ^ 73a [Y2 (gi a]i T3[F2 <8) a]i 

i3iln(2p7^) 1.633 1.633 -0.817 -0.817 -0.652 -0.652 

i3iln(l/J^) 8.281 8.281 -2.070 -2.070 -0.661 -0.661 

^^^In(2^?3^) 2.582 2.582 2.582 2.582 0.206 0.206 

i3iSn(3s7^) 0.000 0.000 2.449 -2.449 0.000 0.000 

^^iSn(2d3^) 4.648 -4.648 -1.549 1.549 -0.618 0.618 

i^iSn(2df^) 5.797 -5.797 2.898 -2.898 0.330 -0.330 

^^^Sh{lgl) 12.470 12.470 -2.494 -2.494 -0.711 -0.711 

i33Sb(2dJ) 5.797 5.797 2.898 2.898 0.330 0.330 

i33Sn(2/7) 9.621 -9.621 3.207 -3.207 0.427 -0.427 

^^^Snhhl) 17.160 -17.160 -2.860 2.860 -0.761 0.761 

2 
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TABLE XV. Magnetic moments of light odd-A nuclei. 



s.p. SM exp. 

"B 3.790 2.872 2.689 

i^N -0.263 -0.307 -0.322 

27 Al 4.790 4.207 3.642 

29p 2.790 1.088 1.235 

2.790 1.252 1.132 

=^^01 0.126 0.634 0.752 



TABLE XVL The renormalization of single particle matrix elements in light odd-A nuclei. 







r 




a 




[Y2 (8) a]i 


T3[Y2 (g) a]i 


iiB 


S.p. 


2.582 


2.582 


1.291 


1.291 


0.206 


0.206 




SM 


3.100 


3.100 


0.773 


0.773 


0.309 


0.309 


13n 


s.p 


1.632 


1.632 


-0.408 


-0.408 


-0.651 


-0.651 




SM 


1.657 


1.657 


-0.432 


-0.432 


-0.598 


-0.598 


27 Al 


S.p. 


5.787 


5.787 


1.449 


1.449 


0.330 


0.330 




SM 


6.164 


6.164 


1.080 


1.080 


0.321 


0.321 


29 p 


s.p 


0.000 


0.000 


1.225 


1.225 


0.000 


0.000 




SM 


0.910 


0.910 


0.314 


0.314 


0.246 


0.246 


31 p 


S.p. 


0.000 


0.000 


1.225 


1.225 


0.000 


0.000 




SM 


0.822 


0.822 


0.402 


0.402 


0.190 


0.190 


33Q 


s.p. 


4.648 


4.648 


-0.775 


-0.775 


-0.618 


-0.618 




SM 


4.361 


4.361 


-0.488 


-0.488 


-0.695 


-0.695 
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progress has inspired the calculations reported here. In our work we employ a PNC nucleon- 
nucleon interaction derived from a vr-, p-, and u;-meson exchange model, providing sufficient 
degrees of freedom to describe fully the five independent S — P amplitudes. The single- 
nucleon, exchange current, and nuclear polarization AM contributions are then evaluated 
with this choice of potential. The end result is an analysis of AM constraints that is fully 
consistent with the existing analysis of Al{p + p) and other hadronic tests of PNC. 

Our results show that the weak meson-nucleon couplings favored by nuclear experiments 
are not compatible with the large AM value extracted from the Cs measurement. The Tl 
AM limit also favors a sign disfavored by theory. Our qualitative arguments about the 
effects of correlations missing from the SM calculations suggest that improvements in the 
nuclear structure are likely to lead to smaller values for the predicted Cs AM, exacerbating 
the current discrepancy. 

The nuclear constraints favor a small value for and isoscalar PNC couplings near the 
DDH "best values." This pattern is puzzling, and suggests that strong interactions modify 
the isospin of weak-meson nucleon couplings in a nontrivial way. The Cs AM result now 
has produced a more confusing situation, one where no one solution satisfies all constraints. 
Hopefully new experiments will provide the redundancy needed to resolve the conflict. In 
the next few years results are expected for the spin-rotation of polarized slow neutrons in 
liquid helium |55] and the asymmetry in polarized neutron capture n + p —>■ d + [56|. 

New AM measurements could also help clarify matters. A more accurate Tl AM mea- 
surement could define the sign of this quantity: while the current band includes zero, it 
favors a sign opposite that predicted by theory. New AM measurements in odd-neutron 
nuclei would have great impact, defining a band in the weak meson-nucleon coupling plane 
roughly perpendicular to the Cs and Tl bands. There are proposals for AM measurements 
on Dy, Fr, and Ba^. 

The accuracy of the Cs AM results sets it apart from any other atomic PNC result: 
it has produced a constraint on a weak radiative correction that, when translated into 
meson-nucleon weak couplings, is as accurate as any direct probe of hadronic PNC. Thus 
the challenge of understanding this special measurement should motivate more theoretical 
work. Furthermore, the implications of this measurement are not necessarily limited to the 
issues discussed in this paper. Our understanding of V(e)-A(A^) interactions also affects the 
interpretation of electron- nucleus scattering experiments like SAMPLE |]5^ , where a similar 
discrepancy between theory and experiment exists and where theoretical predictions also 
depend on a proper treatment of the hadronic weak interaction. Unravelling the puzzles 
presented by these measurements constitutes an important challenge to both theory and 
experiment. 



ACKNOWLEDGMENTS 

This work was supported in part by the US Department of Energy and by the National 
Science Foundation. WH thanks the Guggenheim and Miller Foundations for support during 
completion of this work, and UC Berkeley for its hospitality. 



41 



APPENDIX I: TWO-BODY EXCHANGE CHARGE AND CURRENT 
OPERATORS IN MOMENTUM SPACE 

The total Lagrangian density we are considering is 

L = Cpree + ^PC + ^PNC + ^EM 5 (l) 

with 

Cpree = n\i(^- M)N + \{d,n) ■ (d^n) - ^mln' - ^F^^^ ■ F'^^'^^'^ + \mlp, ■ (T 

- + \mlu,uj^\ (2) 
C^pc = W-kNnN 75r ■ TxN - QpnnN (7^ - i^a^^q'')T ■ p^N 

- 9.NNN\j,-i^a,,qnco^N, (3) 

CpNC = -^N'{rX7r),N + N\hlr-fr + h^, + ^(SraP^ - r ■ p-^))7^75iV 

+ n'Ku^ + hlnu>^)^,^,N + 0(/iJ'), (4) 



e vr X 



Note that we use the Bjorken-Drell [5S] metric exclusively and the DDH definition of weak 
couplings. In these expressions tt, p^, u'^, denote the pion, rho meson, omega meson, and 
photon fields; Fj^^'^''^^ is the field tensor for the designated field; and and are the 4- 

momenta carried by the outgoing meson and photon. F^*' 2'^"* denotes the isoscalar or isovector 

EM form factors, with ^(0) = fI^^ = 1, F^^\o) = p, = -0.12; ^^^^(0) = p„ = 3.70. 
After applying the procedure described in Sec. ^, we obtain the following results. 



A. Pair Currents 



Pair current diagrams are generated by tt, p, or uj exchange, and the nucleon coupling to 
the photon has either a PC or PNC meson-nucleon coupling. Thus there are six cases. For 
charge densities to 0(1/M^) we obtain 

p'!^-PNc = 0, (6b) 
P'-ph = [(1 + Ps){hlm . f(2) + hlT,{l) + i(3r3(l)r3(2) - f(l) ■ r(2))) 
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+ (l + /x.)(/^jr3(2) + /ij + -^r3(2 

H k ■ a(l) X -P2) 

mi 



k ■ a{l) X a(2) 



+ + --^)(r(l)xr(2))3 



X 



/c • ct(2) + — k ■ - P2) 



mi 



" 2v^- 

{P2 -P2y + mi 



+ (1^2), 



p pair 
Pj-PNC 



+ (l--2), 



h 



(1 + - -^)(f(l) x f(2))3a(l) • k-^ 



(27r)35(3)(. 



2^6 



(P2 -nf + ml 



pV-Tc = ^^(^° + ^L^3(2))[(l + + (1 + /..)r3(l)] 



A; • <t(1) X a(2) 
(27r)=^5(3)(...) 



+ (1^2), 



LO pair _ r. 
P-y-PNC — ^■ 



For current densities to 0(1/M^) we obtain 



■7r pair 



2V2M 



f(l) . .12) - .3(l)r3(2))^(l)^|^)Ut_L + (1 „ 2), 



•TT pair /-> 

Jl-PNC — ^1 



■ p pair 
J-1-PC 



egpNN j 
2M \ 



hl{f{l) ■ f(2) + r3(2)) + h\{l + T,{1)) + i(3r3(l)r3(2) 



f(l)-f(2) + 2r3(2)) 



2V6 

^^(2) + 5 (P2 -P2) 



mi 



+ 



2^6 
^(2) • (P2 - P2) 



)(f(l)xr(2)), 



a(l) X -P2) 



a(l) X a(2) 

(27r)3<5(3)(. 



(P2 -P2)^ + m] 



-^ppair _ —^9pNN 
Jj-PNC — 



+ (1^2), 



(6c) 

(6d) 

(6e) 
(6f) 

(7a) 
(7b) 



(7c) 



h^iril) ■ f(2) + r3(2)) + hl(n(2) + rs(lM2)) + ^{3^(1^2) 



2M 

f(l).f(2) + 2r3(2)) 



{P2 - P2Y + ml 



pair ^QwNN / , o 

Jj-PC — 



■u) pair 
Jl-PNC 



2jvf (A!, + AiT3(2))(l+T3(l)) 

-'=''"™(ft°+fti)(l + r3(l))?(l)#>'*'"<- 



(^2 - ^2)^ + ml 



+ (1^2), 



2M 



(P2 -P2)^ + ml 



+ (1^2). 



(7d) 

(7e) 
(7f) 
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B. Transition Currents 



The transition currents can have a 7r7r7, pp7, p7r7, or ujt:^ vertex. In the last two cases, 
the heavier mesons p and u can have either a PC or PNC couphng. Thus there are six 
possibihties. For charge densities to 0(1/M^) we obtain 



= ^?^^[r3(l)r3(2) - f(l) • t{2)]{E; - E, - (E; - E,))a(l) ■ {p[ - n) 



X 



(27r)35(3)(. 



3 



+ (1^2), 



(8a) 



p'"'^ = ^e^,,v^^(/.°-^)(r(l)xr(2)) 
^(2) • {V2 - V-i) 



+ 



-(P2 -P2) ■ (Pl -Pl) 



(27r)35(3)(...) 



+ (1^2), 



(8b) 



P^PC = T/?ff,^^^ (T(l) X f(2))3 \{p[ - p,) ■ {p; - p,) X {p[ + p,) 



2V2M- 



+ «(1 + l^v){^{'^) ■ (Pl - Pl)(P2' - P2) • (Pl - Pi) - ^(1) • (P2' - P2)(Pl - Plf) 

(27r)35(3)(---) 



X 



((pi -pi)^+"^p)((P2 -P2)^+^^; 



+ (1^2), 



pT-pnc = '-^^fi^(Kni) . m + hlr,{2) + ^{3r,{l)M2) - f(l) ■ f(2))) 



(8c) 



X ^(2) • (Ps - P2)o'(l) ■ (Pi' - Pl) X {P2 - P2) 

(27r)3<5(3)(...) 



X 



((Pi -Pi)^ + ^p)((P2 -P2)^ -m2) 



+ (1^2), 



Pu,-PC = 0> 



Poj-PNC — 



2 Mm, 



(^^3(2) + /iir3(l)r3(2))a(2) • (p, -p^) 

(27r)3<5(3)(...) 



X <?(!) • (Pi -Pi) X (p2 -P2)77^ 



((Pi -Pi)^ + ^S)((P2 -P2y + ml) 
For current densities to 0(1/M2) we obtain 

J'^'^" = ^|^[r3(l)r3(2) - r(l) ■ f(2)]((p-' - p^) - -p'i))a(l) ■ {pi - p,) 



(8d) 
(8e) 

+ (1^2). (8f) 



X 77^ 



(27r)35(3)(. 



((Pi - Pi)^ + ml){{p2 - P2Y + ml] 



+ (l--2), 



(9a) 



2M 



- ^) m X f(2))3| [(p-; - Pl) - (pt - p-'2)] (a(2) 



+ 



^(2) • (pj - P2) , ^' 



mi 



(P2 - P2)) • ((P2' +P2) - (Pi' +Pi) - «(1 + /^^;)^(1) X (pi - Pl)) 
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+ 



+ 



ip; + p^) + t{l + X - p^)] {p^ - p,) ■ (a(2) + ^^^^^tS^^JlL{p^ 



ai2) + '-^^^^4^ip;-p.) 



mi 



+ X {p[ -pi) 



{P2 - n) ■ {{P2 + P2) - {Pi + Pi) 
(27r)3<5(3)(...) 



((Pi - Pif + ml){{p2 - P2Y + ml) 



+ (1^2), 



f;jpc - '-^f^0^iril) X f(2))3(p; -pO X (p-' -P2) 



X 



(27r)3(5(3)(. 



Jp-PNC 



{{pi - Pi)^ + ml){{p2 - P2)^ + ml) 



+ (1^2), 



X a{2) ■ - P2)]{ [(Ps' - Pi){K - Pi) - (Pi' - fi){P2 - P2)] X a{l) 
- • (Pi + Pi)(Pi - Pi) X {P2 - P2)| 

+ (1^2), 



(27r)35(3)(. 



((Pi -Pi)^ + ml){{p2 -P2)2 + m2) 



JUJ-PC — ^1 



Juj-PNC 



^^^(/.^3(2) + /iLr3(l)r3(2))a(2) • (p1 - P2) { [(p','^ - p1) (p',' - Pi) 
(Pi^ - P?)(P2' - P2)l X B{1) - a{l) ■ {p[ + pi){p[ - pi) X (P2' - P2) 



X 



(27r)35(3)(. 



((Pi - Pi)^ + ml){{p2 - P2)2 - m^) 



+ (1^2). 



P2 



(9b) 



(9c) 



(9d) 
(9e) 



(9f) 



APPENDIX II: TWO-BODY EXCHANGE CURRENT OPERATORS IN 
POSITION SPACE TO ORDER OF 1/M 



Only the three-current operators are needed for the AM calculation. We keep terms 
to 0{1/M). The following results follow from Fourier transformations of selected terms in 
Appendix A. 



h-ph = • ^(2) - r3(l)r3(2))a(l)5(^)(f - fO- 



Sv^ttM 



+ (1^2), 



(1) 



pair (^9pNN j 



hl{f{l) ■ f(2) + r3(2)) + hl{l + r,{l)) + ^{3r,{l)r,{2) - f(l) • f(2) 



+ 2r3(2))] (a(2) + x/2^(l + + (;;^)r2(f^.) ® a{2)]i) 



+ {K - X r(2))3(^(l) X a{2) - ^{1 + 



3 3 , 
+ 7 -) 



{nipry 
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X [F2(^^r) ® X a(2)]i+iv^(l + -^ + -^)[F2(l^r-)® [a(l)®a(2)]2]i' 



■p pair 



hi 



hl{T{l) ■ f(2) + r3(2)) + /.J(r3(2) + r3(l)r3(2)) + ^(3r3(l)r3(2) 



-^ui pair _ ^9u)NN /j^Q 
Jj-PC 



•u) pair 
J-y-PNC 



- r(l) ■ f(2) + 2Ts{2))]a{l)5^'\x - Xr)- + (1 ^ 2), 

X [y2(^^r-) ® <?(2)]i)(5(^)(f - fi)— — + (1 ^ 2), 

_ —eguNNf,o , ,1 



AttM 



{hi + hl){l + T,{l))a{l)5^'\x-x,)— + (1 - 2), 



J.., = -|^^^(f(l) . f(2) - r3(l)r3(2))a(l) • Vi(Vi - V2) 



16V27rM 



(2vr) 



ik-(x—R) 



^ ^^9pNN ,,Q _ 



-1/2 

)(r(l)xr(2))3 



+ (1^2), 

a(2) • V2 



(2) 
(3) 

(4) 
(5) 

(6) 



(a(2) - 



V2) • (Vi - V2 



+ i{l + pi,)a{l) X Vi)(Vi - V2) - Vi • (a(2) - ^i^LZ^V2)(Vi 



+ i(l + X Vi) + V2 • (Vi - V2 + i(l + x Vi)(a(2) - ^l^LZ^Vs) 



X 



/ 



-Lor 



^ iH^-ii) f"^ dae^ak-r"—- + (1 ^ 2), 
{iTiy^ J-i/2 Lp 



^ 8v27rmo ^ (27r)'* 7-1/2 



+ (1^2), 



(7) 



(8) 



where R = {xi + X2)/2, f = Xi - X2, r = \f\, L^(p) = (m^(^) + P{1 - a^)y/^, L 



^m^+nv _|_ Q,^^2 _ ^2^ _|^ ^2|-i _ ^2^^1/2^ Operation of V should be understood as 

V F{- ■ ■) = F{- ■ ■)V — VF(...) . For transition currents, the full Fourier transformation is 
not easily evaluated in the general case, so we leave the integration undone. 



APPENDIX III: AN EXAMPLE OF FERMI GAS ONE-BODY AVERAGES 

Here we describe how effective operators are obtained in the Fermi gas model by one- 
body averages. Most of the discussion is general, though we use the simple tt pair current 
when specific examples are needed. 
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TABLE XVII. One-body averaged spin operators. 



2-body 1-body direct 1-body exchange 
i 2 1 
(?(!)+ (7(2) 2a 2a 
a{l) - a{2) 2a 
a{l) X a{2) 2ia 



TABLE XVIII. One-body averaged isospin operators. 



2-body 1-body direct 1-body exchange 

1 {Op + 0n) {9p + 9n)l + iep-9n)f 

r(i)-n2) (ep-OnH {0p + 0n)l-iep-9n)^ 

r3(l)r3(2) i9p-enH {Op + 9^)1 + {9p - 9n)^ 

r3(l) + r3(2) {6 

p 9^1 ) + {9p + 6n)T'i (dp - 9n) + {9p + 6n)T3 

r3(l)-r3(2) _(^p_^„) + (0p + 0^)^3 

(f(l) X f(2))3 -ii9p - On) + i{9p + 0„)r3 

3T3(l)r3(2)-f(l)-f(2) 2(0p-0„)r3 2(0p - e„)r3 



An effective one-body operator is obtained by performing a mean-field-like sum over the 
direct and exchange terms 

< a|0(^)|/5 >=Y.< "710(^^1/57 > - < a7|0(2)|7/5 >, (1) 

7 

where the sum extends over occupied core states. In the nonrelativistic Fermi gas each s.p. 
state is a direct product of space, spin, and isospin components 

\a) = \p{a)) ® \^ms{a)) ® \^mt{a)). (2) 

Thus the wave function factors, allowing space, spin, and isospin sum to be performed 
independently. 

The spin and isospin averages for common operators are easily done. The results are 
displayed in Table pOTTll and pCVTT^ . 

Turning to spatial averages, we first consider pair currents. The spatial parts of these 
operators take one of the generic forms: i) f{r)5^^\x — Xi) or ii) f(r)5^^\x — X2), where /(r) 
is a function of r = |r| = |xi — a;2|. Therefore, the direct average is 

E(P-P7l/W^^'H^-^i)b">,P7) = e-^(P^'^-^>)-^^E J d^fir), (3) 

P-Y Pi 

Pi Pi 
and the exchange average 
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P7 Pi 



Pi 



(5) 
(6) 



Pi 



Now e-^(P--P'3)-^ gives a 5^'^\x ~ Xi) for a first quantized operator in position space, i.e., 
{Pa\5'^^\x — Xi)\pi3) = e~*'^'°"~^'3)'^. We specialize the remaining integration to the vr pair 
current where /(r) has the Yukawa form e~™""/r. 



r- = — 

TV 

rriTr 



An 



r p + 



2 ■ 



(7) 
(8) 



After performing the sum over the Fermi sphere by using the quasi-continuum hmit, J2p^ 
dpy It dn{py), we find 



El 



Pf P 



Pi 



1 



P 



■W{pa{/3),m^ 



(PaW - P^Y + ml 2ml 
where the W' function represents the full result after the volume integration, 

■I - pa 



(9) 
(10) 



3m? { 

W\p^,m^) = 2-2m, 
4^Pf [ 



+ —[l-pl+ml] In 



arctan ( — z — -) + arctan ( — z — -) 



(11) 



.(1 + m2. 

with all the tilde quantities normalized by the Fermi momentum, i.e. X = X/pp. As shown 
in Fig. [ly, the W function varies very slowly as pa runs from to 1. Therefore, it is 
reasonable to replace this quantity by its average value 

Combining the spatial result with the spin and isospin factors from the tables the yields 
the TT pair current one-body averged form 

h-Tc = ^^((^P + - - 0n)T,)a'^^{W'^-^)6^'\x - xd + (1 - 2) (12) 



(13) 



Other currents are similar, though generally more tedious. 



APPENDIX IV: ONE-BODY FERMI-GAS AVERAGED CURRENT OPERATORS 

We list the relevant one-body Fermi gas averaged operators in momentum space. 
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FIG. 10. The smoothness of W ^'^> as a function of Po 



(1) 



+ 



6Mm2 



(2) 



V 4 y 



hi 



2Mmy " 



+ 



2Mml 



3h 



hi 



— 



'15 1 
4 m? 



CO 



■u) pair 
J-f-PNC 



•TTTT'y 



^^ihl + hl){l + T,)9, p{2(VF'(-))a - y f ® K], ® a],}, (4) 

-^a.NN^^o^ + hl){i + r3)((^„ + e,) pa - e, p{W'^-^)a), (5) 

m + ^p) + {On - Op)Ts) p{{Y,){a ■ k)k + {Y^){a ■ K)K + {Y^){pl)a}, (6) 



2Mml 

—(^g-KNNf-K 



)((^n - Op) + {9n + ^p)t3) p[{Il){k'')a + {l2){pl)5 

+ (73)(? • k)k + {h){a ■ K)K - i{h)k xK+ -\({J^){k^)a + {J2){k ■ Kfa 



+ {J^){plP)a + {J,){p%)a + {J,){a ■k){P)k + {J,){a-K){k-K)k 

+ (J7)(p|)(? ■k:)k + {Js){a ■ K){P)K + {J^){py){a ■ K)K + (Jio)(? ■ ^)(^ ■ K)K 



49 



TABLE XIX. Average weighting functions. (Note: First number refers to proton part and the 
the second the neutron part) 





\ / 


\ / 


\ 1 


\ -L / 






0.30/0.26 


0.90/0.88 


0.99/1.02 


0.19/0.23 


0.0039/0.0033 


0.0173/0.0111 


0.0144/0.0103 
















2.25/2.18 


10.09/9.58 


-1.53/ - 1.46 


21.30/19.61 


10.18/9.56 


0.42/0.39 


1.10/1.02 
















4.11/3.91 


1.43/1.32 


-0.83/ - 0.77 


-2.68/ - 2.48 


2.57/2.40 


4.25/3.89 


-8.21/ - 7.56 


(^10) 


(^11) 


(^12) 










2.53/2.32 


3.26/3.03 


-11.16/ - 10.42 


1.43/1.15 









- «(Jii)(P)A; X - t{Ji2){pl)k X K 

Qp(n) is a projection operator: when Qp(n) acts on {p (nuclear density), "pp (Fermi momen- 
tum), (X) (averaged weighting function such as (VT '•'^^))}, the results should read {pp(n), 
(Xp(„))}. The conversion rules from momentum space to position space are simple, 

1 _^5(3)(^_^^)^ (9) 
^ = P/3 - Pa — ^ -^[V, - (10) 

K = _,{v„5(3)(^ _^^)}^^^. (11) 



For ^^^Cs, p^^"^ ~ 260(300) MeV. The average weighting factors are given in Table XIX . 
For ^°^T1, these numbers are almost the same. We do not distinguish between the p and u 
masses in evaluating weighting functions. 
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